Novel Phases and Reentrant Melting of Two Dimensional Colloidal Crystals 
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We investigate two-dimensional (2d) melting in the presence of a one-dimensional (Id) periodic 
potential as, for example, realized in recent experiments on 2d colloids subjected to two interfering 
laser beams. The topology of the phase diagram is found to depend primarily on two factors: the 
relative orientation of the 2d crystal and the periodic potential troughs, which select a set of Bragg 
planes running parallel to the troughs, and the commensurability ratio p = a' /d of the spacing a 
between these Bragg planes to the period d of the periodic potential. The complexity of the phase 
diagram increases with the magnitude of the commensurabilty ratio p. Rich phase diagram, with 
"modulated liquid", "floating" and "locked floating" solid and smectic phases are found. Phase 
transitions between these phases fall into two broad universality classes, roughening and melting, 
driven by the proliferation of discommensuration walls and dislocations, respectively. We discuss 
correlation functions and the static structure factor in these phases and make detailed predictions of 
the universal features close to the phase boundaries. We predict that for charged systems with highly 
screened short-range interactions these melting transitions are generically reentrant as a function of 
the strength of the periodic potential, prediction that is in accord with recent 2d colloid experiments. 
Implications of our results for future experiments are also discussed. 
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I. INTRODUCTION 
A. Motivation and background 

Two-dimensional (2d) melting and the mathemati- 
cally related systems, such as for example the normal-to- 
superfluid and planar paramagnet-to-ferromagnet transi- 
tions in films (described by the 2d XY model) are striking 
examples of the increased importance of thermal fluc- 
tuations in low dimensional systems. In contrast 
to their bulk, three- (and higher-) dimensional analogs, 
where, typically, fluctuations only lead to quantitative 
modifications of mean-field predictions (e.g., change val- 
ues of critical exponents), here the effects are qualitative 
and drastic. Located exactly at the lower-critical dimen- 
sion (di c = 2), below which the distinction between the 
high and low temperature phases is erased by fluctua- 
tions, two-dimensional melting can proceed via a sub- 
tle, two-stage, continuous transition, driven by unbind- 
ing of topological defects (dislocations and disclinations). 
This mechanism, made possible by strong thermal fluctu- 
ations, therefore provides an alternative route to a direct 
first-order melting, argued by Landau's mean-field anal- 
ysis H to be the exclusive scenario. 

Despite of its long history, dating back to the work of 
Kosterlitz and Thouless || , Halperin and Nelson, H , and 
Young H (KTHNY) (which in turn built on a large body 
of ideas dating back to Landau and Peierls 0), interest 
in 2d melting and related problems persist. On the theo- 
retical side this, in part, is due to the fact that the theory 
of 2d melting is an unusual example of a nontrivial and 
quite exotic critical point that lends itself to an asymp- 
totically exact description. Furthermore, the KTHNY 
class of transitions (2d melting and related disordering 



of a 2d XY model) provides a rare example of a ther- 
modynamically sharp phase transitions between phases, 
both of which lack long range order. || 

Although evidence for defect driven phase transitions 
has appeared in numerous experiments on liquid crys- 
tals H and Langmuir-Blodgett films [nj, finding simple 
model systems which exhibit these phenomena in exper- 
iments or simulations has proven to be more controver- 
sial. Some system parameters appear to fall in the range 
in which instead it is discontinuous melting that converts 
a solid directly into a liquid. However, it appears, that 
two-stage continuous melting has been recently experi- 
mentally observed by Murray etal. |pd| and Zahn etal. 
p"2[ in beautiful melting experiments on two-dimensional 
colloids confined between smooth glass plates and super- 
paramagnetic colloidal systems, respectively. In these 
experiments, an orientationally quasi-long-range ordered 
but translationally disordered hexatic phase jBj was ob- 
served. This phase, intermediate but thermodynamically 
distinct from the 2d solid and isotropic liquid, is an im- 
portant signature of defect driven two-stage melting. In 
these two-dimensional colloids, particle positions and the 
associated topological defects can be directly imaged via 
digital video-microscopy, allowing precise quantitative 
tests of the theory. Colloids are thus ideal experimental 
model systems to explore the details of two-dimensional 
melting and related phenomena, many of which are the 
focus of the theory presented here. fl3|| 

Soon after the initial development of the theory of 
two-dimensional melting, theoretical efforts turned to the 
studies of the effects of substrate, an important ingrc- 
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dient in many physical systems. These studies J14 
uncovered a rich phenomenology stemming from the in- 
terplay between the underlying periodic substrate and 
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a quasi-long-range ordered solid film interacting with it. 
While many experiments have been undertaken, with a 
krypton film on a graphite substrate (see e.g. Ref. Q for 
a review) being one of the best studied, these systems are 
far from ideal in exploring this rich phenomenology, be- 
cause of the lack of substrate tunability; in these systems 
it is difficult to change the substrate period, dimension- 
ality and pinning strength. 

A series of pioneering experiments by Chowdhury, Ack- 
erson and Clark fl6|| constituted an important new de- 
velopment. In these studies strongly interacting colloidal 
particles, confined to two dimensions, were subjected to 
a one-dimensional periodic potential, induced by the in- 
terference fringes from two laser beams crossed in the 
sample. The light-induced polarization in these micron- 
size dielectric particles interacts with the laser electric 
field, leading to a radiation pressure force, jL7| directed 
toward the regions of high laser intensity, the antinodes 
maxima in the laser standing wave pattern. 

One of many interesting phenomena discovered by 
Chowdhury etal. is the fixed-temperature freezing tran- 
sition driven by increasing the strength of the laser po- 
tential, dubbed "light induced freezing" (LIF). Qualita- 
tively, LIF is due to the suppression of thermal fluctua- 
tions transverse to the imposed periodic pinning laser po- 
tential. This intuition is also supported by a more quan- 
titative analysis based on Landau's free energy expansion 
in the translational order parameters (density Fourier 
modes) pc i: with {G^} the three smallest reciprocal lat- 
tice vectors of a triangular lattice. In the simplest geome- 
try, with e.g., Gi commensurate with the laser potential, 
(pd) is trivially induced by the potential even in the liq- 
uid phase. Such finite (pd) then converts the Landau's 
cubic coupling pgi/Og 2 Pg 3 , (which, in mean-field theory, 
is responsible for melting always being first-order) into a 
simple upward shift in the melting temperature for the 
only remaining critical mode ip = pc 2 — pg 3 . Not sur- 
prisingly the resulting Landau expansion contains only 
even powers of this complex order parameter ip, which 
therefore generically orders via a continuous transition 
in the XY universality class. Hence, within the mean- 
field description discussed by Chowdhury etal. [|l6|, one 
expects to reach a tricritical point upon increasing the 
light intensity, beyond which the LIF transition becomes 
continuous. 

However, because of the dominant role of thermal fluc- 
tuations in two-dimensional systems, such "soft-spin" 
Landau expansions in order parameter amplitudes (and 
the related density functional theories |l^]) will have 
difficulties to capture the subtleties of the continu- 
ous topological phase transitions possible in these two- 
dimensional systems. Unfortunately, results from Monte- 
Carlo simulations are inconclusive. Although earlier sim- 
ulations fl9[| claimed to have found a tricritical point 
at intermediate laser intensities, consistent with density 
functional theory, recent studies from the same labora- 
tory [^0| refute these results. These difficulties are per- 
haps unsurprising, given that even much larger scale sim- 



ulations have, so far, failed to completely resolve the na- 
ture of 2d melting, even without a periodic external po- 
tential. |nj 

An alternative (but complementary and in principle 
equivalent) "hard-spin" defect description (with order 
parameter amplitude fluctuations represented by defect 
cores), extended to include a one-dimensional periodic 
pinning potential may be necessary to correctly cap- 
ture the rich phenomenology of the early experiments by 
Chowdhury etal. |l6| and the recent ones by Wei etal. |2^] 
and others. Developing such a theoretical framework 
and exploring its details to interpret these experiments 
is the goal of the work presented here. 

Our interest in this problem was stimulated by the 
experiments of Wei etal. p2fl , which extended the light- 
induced melting experiments to higher laser intensities 
than those studied in Ref. One other notable dif- 

ference is that in contrast to the strong long-range in- 
teraction of unscreened charged colloids in highly deion- 
ized solution, |l6| in the Wei etal. experiments col- 
loidal particles were interacting via a short-ranged Debye 
potential, with ions in the solution screening the long- 
ranged Coulomb interaction. In addition to the light- 
induced freezingjObserved at low light intensities, the 
authors of Ref. [|22f discovered a reentrant melting phe- 
nomena, "light-induced melting" (LIM), driven by the 
increased strength of the laser-induced one-dimensional 
periodic potential. As discussed below, this fascinating 
reentrance phenomena generically emerges from our the- 
oretical analysis in the limit of a short Debye screening 
length. 

The goal of this paper is to investigate the phenom- 
ena of two-dimensional melting in the presence of a one- 
dimensional periodic potential, and to answer many basic 
questions stimulated by these recent experiments. What 
is the nature of such melting transition, if not preempted 
(as it can always be) by the first-order transition? More 
generally, how is the standard phase diagram for 2d melt- 
ing on a homogeneous substrate (which includes the 2d 
crystal, hexatic and liquid phases) modified by the pe- 
riodic laser potential? Which of the phases survive the 
light field and what new ones emerge in its presence? 
The answers to these and many other questions, pro- 
vided below, lead to results consistent with experimental 
observations, and have many testable consequences for 
the possible future experiments. 



B. Summary of the results 

Even in the liquid phase at high temperatures the laser 
interference fringes, which we choose to run along the x- 
axis, induce a periodic density modulation in the colloidal 
liquid. As a consequence the static structure function 
5(q) displays Bragg peaks at K„ = n{2i:/d)y, the in- 
teger multiples (n E Z) of the reciprocal lattice vector 
K = (27r/d)y of the imposed one-dimensional periodic 
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potential with trough spacing d. j24| The liquid phase 
density exhibits a finite linear response to such a periodic 
perturbation with amplitude Uk, which is proportional 
to the input laser intensity J; n . This is consistent with 
the observations of Chowdhury etal. 16 who found the 
scattered laser intensity, 7 ou t, at these directly induced 
Bra gg p eaks, to scale as a cube of the input laser power 
/i n . Pq] These explicitly induced features of the modu- 
lated liquid persist throughout the phase diagram, with 
the additional structure emerging as a result of numer- 
ous spontaneous symmetry breakings, which we discuss 
below. 

The laser-induced periodic potential also explicitly 
breaks continuous 2d rotational symmetry down to Z2 
symmetry (rotations by n). Consequently, the one- 
dimensional periodic potential induces nematic, square, 
hexatic and higher orientational harmonics long-range or- 
ders, respectively characterized by a 2n-atic bond ori- 
entational order parameter ip2n — (e t2n0 ^), which, in- 
dependent of any other details, are nonzero throughout 
the phase diagram. Therefore, in particular, the laser 
potential eliminates the continuous transition from an 
isotropic liquid to a hexatic liquid phase, expected in 
two-dimensional liquids in the absence of an external po- 
tential. ||] This situation is analogous to a ferromagnet 
in a magnetic field, where the qualitative distinction be- 
tween paramagnetic and ferromagnetic phases is erased 
by the external magnetic field, with both phases display- 
ing a finite induced magnetization. 

Since the hexatic orientational order is explicitly in- 
duced by the laser potential, it must vanish as the laser 
field is tuned to zero. Analogously to a power-law vanish- 
ing of the magnetization with an external magnetic field 
in a ferromagnet at its critical point, we predict that at 
low input light intensities, Ji n , the orientational order pa- 
rameter vanishes as a universal power of ii n , 



with \/5^ 6 = 6 Pqj in the liquid phase and 
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in the hexatic phase, where fj 6 is the exponent describ- 
ing the algebraic decay of bond orientational order in the 
absense of the laser- induced periodic potential. || We ex- 
pect iI>q to approach a nonzero /; n -independent constant 
in the solid phase, consistent with the spontaneous long- 
range hexatic order of the 2d crystal, even in the absence 
of a periodic potential. 

All other details of the phase diagram and the proper- 
ties of the phases for our system strongly depend on the 
level of commensurability between the two-dimensional 
colloidal crystal, in the absence of the laser field, and the 
one-dimensional periodic potential that it induces. This 
in turn is determined by two ingredients: (i) the orien- 
tation of the triangular colloidal lattice relative to that 



of the periodic potential troughs, which selects a set of 
Bragg planes that run parallel to the troughs, (ii) the 
commensurability ratio of the spacing a' between these 
Bragg planes to the period d of the laser potential, de- 
fined by p = a! I d. In this paper we will primarily focus 
on the commensurate case defined by p € Z and defer 
the rich phenomenology of the incommensurate case and 
the commensurate-incommensurate transitions to a later 
publication. |27| 

For these commensurate densities, independent of the 
order of commensurability, p, at the lowest temperatures 
we always find that our system freezes into a novel type 
of a crystal, which we call a "locked floating solid" (LFS). 
This phase derives its apparently contradictory name 
from its novel highly anisotropic properties: while the 
solid is pinned transversely to the troughs of the peri- 
odic potential, executing only massive optical phonon- 
like excitations in that direction, it is able to slide freely 
along the potential minima with acoustic phonon excita- 
tions within the troughs. Upon integrating out the mas- 
sive ity-modes and using standard renormalization group 
methods [^) to eliminate bound dislocation pairs in the 
LFS phase, we are left with a free energy with temper- 
ature and potential strength dependent effective elastic 
constants, 



- I d 2 r{K cS (d x u x ) 2 + fi cS (d y u x ) 2 } . (1.3) 



The structure function of LFS is quite unusual. Like the 
high temperature modulated liquid discussed above, the 
LFS displays a set of delta-function Bragg peaks (reduced 
by the Debye- Waller factor) located at the multiples of 
the laser potential reciprocal lattice vector K = (27r/d)y, 
which coexist with other spontaneously induced Bragg 
and quasi-Bragg peaks. 

The more detailed properties of the LFS and other 
phases exhibited by our system, strongly depend on the 
choice of the infinite set of colloidal crystal orientations 
relative to the light interference fringes. While we will 
explore these numerous possibilities in their full general- 
ity in the main body of the manuscript, in this subsection 
we summarize our results only for the simplest orienta- 
tion studied in the experiments of Refs. ]lq^2]| , in which 
the periodic potential troughs run parallel to the primary 
Bragg planes. |||] 

Experimentally, we expect our system to display a con- 
siderable amount of irreversibility, with the choice of the 
relative orientation highly dependent on the way the sys- 
tem is taken into the crystal state: if the laser potential 
is turned on in the liquid phase (field-cooled), the crys- 
tal will freeze into the lowest energy orientation consis- 
tent with the imposed colloidal density (or the chemical 
potential) and laser fringe spacing; |29| in contrast, in 
zero-laser-field cooling experiments, an already formed 
crystal may be unable to reorient significantly, and will 
therefore lock into a metastable orientation, determined 
by the plane of the two interfering laser beams. 
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FIG. 1. Triangular lattice with lattice constant a subject 
to a periodic potential (maxima indicated by dashed lines) 
for pd = a' with a' = \/ia/2 and p — 2. Also shown is the 
low energy dislocation with Burgers vector b parallel to the 
corrugation of the external potential. 

Once we focus on the primary orientation, illustrated 
for p — 2 in Fig.[l], the phenomenology of our system is 
completely determined by the integer commensurability 
ratio p. As we will show, for commensurate densities, 
our system admits three phase diagram topologies, cor- 
responding to the three ranges of the values of p: (i) 
p = 1, p0| (ii) 1 < p < p c , and (iii) p > p c , with the 
critical value of p c ~ 3.7 for the primary orientation. 



1. Commensurability ratio p = 1: 



isotropic liquid phase is absent in the presence of a peri- 
odic potential, this phase diagram contains only two ther- 
modynamically distinct phases at finite Uk- the mod- 
ulated liquid and the simplest p = 1 "locked floating 
solid" (LFS). We can estimate the order of magnitude of 
the transition temperature between the LFS and liquid 
phases in terms of micros cop ic elastic constants (similar 
to those appearing in Eq |l.3| ) as follows: In the limit of 
strong laser potential the particles are confined to a par- 
allel array of equally spaced Id channels of spacing d, 
illustrated in Fig||. If u n {x) is the particle displacement 
field along the n-th channel, we can write the energy of 
these weakly coupled one-dimensional rows of particles 
as 



/ a \2 
P" — ; cos 



— (u n+ i(x) - u n {x)) 
a 



(1.4) 



where K is the bare compressional elastic modulus within 
each channel and /i is the microscopic coupling between 
the channels determining the shear modulus of the 2d 
system. |3l| l At high temperatures or weak microscopic 
coupling /i, the colloid decomposes into an orientation- 
ally ordered two-dimensional liquid of decoupled one- 
dimensional channels. At temperature T, the phonon 
fluctuations within a channel then grow according to 



For p = 1, we find the phase behavior of the 2d 
colloidal system as summarized by the phase diagram il- 
lustrated in Fig.|[ 




FIG. 2. Schematic phase diagram for a primary commen- 
surate orientation with commensurability ratio p = 1. Th 
indicates the transition temperature from the hexatic to 
the isotropic liquid phase at Uk = 0. Insets: Schematic 
structure functions in the various phases. The x's indicate 
delta-function Bragg peaks and the shaded circles algebraic 
peaks. 

Because the sharp distinction between the hexatic and 



(\u n (x) - u n (0)\ 2 ) = 



(1.5) 



as can be seen from the equipartition theorem. Upon 
choosing x such that the root mean square phonon fluc- 
tuations equal the intrachannel particle spacing a, we de- 
termine a translational correlation length ^t(T), which 
diverges at low temperatures 



K 



Kd 



(1.6) 



n=3 
n=2 
n=l 



(I 
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FIG. 3. Colloidal particles in channels, labeled by n, 
with intrachannel compressional modulus K and interchan- 
nel shear coupling fi. 
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The channels will couple to form a coherent 
two-dimensional LFS when the effective coupling 
di;T(T)[i(a/2iTd) 2 between correlated Id regions of size 
£t(T) surpasses the thermal energy ksT which decorre- 
lates the Id channels. We associate this characteristic 
temperature with the melting temperature T m of LFS, 
which is therefore given by 



knT„ 



const, a 



(1.7) 



A similar argument leads to the estimate for freezing 
into the three-dimensional locked floating solid phases 
discussed by Carraro ]3^] for rare gas atoms adsorbed 
into bundles of carbon nanotubes. As we describe in Ap- 
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pendix ^], in terms of the weakly coupled model, Eq, 
freezing into LFS takes place at strong coupling /i, and 
therefore does not allow a rigorous renormalization group 
treatment of the transition. Nevertheless an approximate 
variational treatment is possible and is presented in Ap- 
pendix |A|. 

Instead, here we take an alternative route to the study 
of the LFS melting and other tra nsiti ons by working 
within a continuum clastic model Eq.1.3, which is equiv- 
alent to the strong coupling ( betw een the channels) limit 
of the descrete model in Eq. 1.4. Such approach allows 
a more refined and asymptotically exact renormalization 
group analysis (presented below), within which we find 
that for p = 1 the melting of the LFS phase is in the uni- 
versality class of the XY model, and is driven by unbind- 
ing of dislocation pairs with Burgers vectors b = ax along 
the troughs of the periodic potential. Consequently, in 
contrast to the conventional 2d melting transition, at the 
melting temperature T m , we predict a universal ratio of 
the jump in the geometric mean of the long wavelength 
effective shear and bulk moduli, ii e g(T~) and K^T^) 
(describing the elasticity of the LFS phase) to T m , 



K e s(T m ) fi c s{T n 



8tt 

TbP ■ 



(1.8) 



This is in agreement, up to constants of order 1, with 
the rough estimate of the melting temperature, Eq.1.7 
sketched above, and with the variational method pre- 
sented in Appendix^. The most striking feature of the 
p = 1-LFS melting transition is the shape of the phase 
boundary T m {Uk), whose universal features guarantee a 
generically reentrant melting, under conditions such as 
the experiments of Wei etal. |22|. At low light intensi- 
ties, i.e., small Uk, we find that the melting curve has a 
universal, cusp shape: 



T m {U K ) ~ T m (0) + [Hk B T m /U K )} 



-i/v 



(1.9) 



with v w 0.36963. On the other hand for large Uk , i.e. , 
for ksT m (UK)/UK <C 1, we find that for short-range 
particle interactions {na > 5.8), T m (Ujc) generically in- 
creases with decreasing amplitude Uk of the periodic 
modulation, according to 



T m (U K ) = T~ 1 



5[(«a) 2 -31] 
64?r2 
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(1.10) 



thus implying reentrant melting for a band of tempera- 
tur es as a function of potential strength (see Fig. ||). In 
Eq.1.10 above, k is the inverse of the Debye screening 
length, tunable by adjusting the solution salt concentra- 
tion, and = T m (UK ~ > oo), which, for the system 
studied in Ref. ]22]| , we estimate to be approximately 
1.3T m (L^ = 0). 

The structure function for the p = 1-LFS, illustrated 
in Fig.^ is also quite unusual. In addition to the set of 
Bragg peaks, directly induced by the laser field, S(q) also 
displays an independent set of quasi-Hv&sg peaks at the 
off-tfy-axis reciprocal lattice vectors G, || 



S(q) 



1 



q-G 



1 2— TIG 



(1.11) 



which distinguishes the LFS from the modulated liq- 
uid state. The corresponding density-density correlation 
function Cg (r) = (/? G (r) J o G (0)) for reciprocal lattice vec- 
tors with G x ^ shows a power-law decay 



Cb(r) 
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(1.12) 



where (i c g an d K c k are the effective shear and bulk elastic 
moduli in Eq.1.3 for the deformations along the troughs 
(cc-axis) of the periodic potential. The exponent r/c de- 
pends on the relative orientation of the colloidal crystal 
and the troughs. Unlike conventional 2d melting 0, it is 
universal at the melting transition, and is given by 



V* G = V G (T~) = (G • b/47r) 



(1.13) 



where b is the smallest allowed Burgers vector in the 
trough direction. For the primary orientation, illustrated 
in Fig.[j], with b = a, the exponent characterizing the al- 
gebraic order in the off-axis peaks (see Fig. |^) closest to 
the Oy-axis is = 1/4; for the next row of peaks with 
G x — 47r I a we find 77^ = 1, consistent with the algebraic 
decay obse rved in Ref. p^] (for a more detailed discussion 
see section VII ). 

Our analysis also makes exact predictions for the struc- 
ture function peak amplitudes in the limit of low laser 
intensity. Simil ar to the hexatic orientational order pa- 
rameter -06 1 Eq.1.1, the translational order parameter, 
defined by Mr„ = (pk„ ) , is induced by the periodic po- 
tential throughout the phase diagram. However, in con- 
trast to the liquid phase, where it vanishes linearly with 
Uk, in the crystal phase for T < T m (0), we find 



Mi 



K„ 



U K 



\1/S M 



(1.14) 



with 5m defined in analogy with the critical exponent at 
a ferromagnetic critical point, 
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1/5 



M 



and where 



k^,T 3/i + A 
4?r fi(2fi + \y 



(1.15) 



(1.16) 



is the exponent with which the real-space density-density 
correlation function decays in a 2d crystal without a 
substrate potential. || We therefore predict that for 
T < T m (0) the intensity of the on-g y -axis Bragg peaks 
vanishes as an exact power of the input laser intensity 



/i n , according to 



iout(K n )~ |<PK„)| 2 / ir 



I: 



1 + 2/&M 



(1.17a) 
(1.17b) 



In contrast, we predict the intensity I ou t of the off-axis 
quasi-Bragg peaks, labeled by a reciprocal wavevector G, 
to vanish as 



/ „t(G) ~ Ii 



l+2f7G/(4-f7 G )^2-(j7 G -i7 G ) 



where 



f) G = rj G (l-G 2 JG 2 



(1.18) 



(1.19) 



and L is the system size. 

We can also define the translational correlation length 
by the widths of the off-g^-axis Lorentzian peaks in the 
structure function. As the melting temperature T m is ap- 
proached from above, given the XY nature of the p = 1- 
LFS melting transition, we expect the correlation lengths 
parallel and perpendicular to the troughs to diverge ac- 
cording to 



,c/\T-T m \ 1/2 



(1.20) 



where c is a constant of order unity. 



2. Intermediate commensurability ratios: 1 < p < p c 

For 1 < p < p c , the phase diagram, illustrated 
in Fig.^, generically includes an additional symmetry- 
allowed "locked smectic" (LSm) phase. 

The LSm distinguishes itself from the modulated liquid 
by spontaneously breaking the liquid's discrete transla- 
tional symmetry by d down to translations by pd. |3^ 1 
In contrast to the LFS, however, the LSm exhibits only 
short-range correlations between colloidal positions lying 
in different troughs, and therefore does not resist shear 
deformations for displacements along the potential min- 
ima. Correspondingly, as illustrated in Fig.||, the struc- 
ture function of the LSm phase displays spontaneously 
induced Bragg peaks at K„/p, in addition to the Bragg 
peaks at K„, directly induced by the laser interference 



fringes. For 1 < p < p c , the LFS also displays these 
spontaneous Bragg peaks on the q y axis at q = ~K n /p. 




FIG. 4. Schematic phase diagram for a primary commen- 
surate orientation with commensurability parameter in the 
range 1 < p < p c (the case p = 2 is shown here). Thin lines 
indicate continuous phase transitions. The thick line between 
the LFS and modulated liquid phase is most likely a first or- 
der phase boundary. Insets: Schematic structure functions. 
As in Fig^, the x 's indicate delta-function Bragg peaks and 
the shaded circles algebraic peaks. 

Symmetry dictates that the freezing of the modulated 
liquid into the LSm is in the p-state clock model uni- 
versality class. Also, similar to the melting of the p = 1- 
LFS, we hnd that the 1 < p < p c -LFS melts into the LSm 
through a transition in the XY universality class and will 
therefore also exhibit the usual Kosterlitz-Thouless phe- 
nomenology. Pi We have also added in the phase diagram 
the possibility of a direct transition from the LFS to the 
modulated liquid at intermediate potential strength. We 
expect this transition to be different than the LFS-liquid 
transition for p — 1. Whereas the p = 1 transition is in 
the XY universality class, for 1 < p < p c the LFS-liquid 
transition is associated with simultaneous loss of XY and 
descrete (Ising for p = 2) order. Because at this latter 
transition two unrelated symmetries are simultaneously 
restored, we expect it to be hrst order. At the multicrit- 
ical point, where the liquid, LFS and LSm phases meet, 
the critical behavior is presumably described by a two- 
dimensional compressible Ising model (for p — 2) [Q of 
the form 



Hi 



I-XY 



[u,S] 



d 2 r 



(VS) 2 



1 



S 2 



vS 4 



+ / d 2 r (j x d x u x + -i y d y u x ) S 2 



Hlfs [u] 
(1.21) 



S is a continuous Ising order parameter, that distin- 
guishes the LSm from the liquid, ^ x ^ y are "magnetoe- 
lastic" parameters, which couple elastic strain to "mag- 



G 



netization" S, and where the parameters of the model 
are tuned to the tricritical point at which both order pa- 
rameters vanish simultaneously. It would be interesting 
to study the properties of such tricritical point, which to 
our knowledge has not been previously explored. 



3. Large commensurability ratios and floating phases: p > p c 

For these higher values of p, the complexity of the 
2d colloidal phase diagram (displayed in Fig.||) further 
increases, allowing two new phases, the "floating solid" 
(FS) and the "floating smectic" (FSm). 

The new phases are distinguished from their "locked" 
counter parts, the LFS and LSm, by their ability to slide 
(float) across the troughs of the periodic potential; tech- 
nically, the periodic potential is irrelevant (in the renor- 
malization group sense) and therefore can be treated per- 
turbatively inside the FS and FSm phases. As illustrated 
in Fig.|5|, all the spontaneously induced structure function 
peaks of these floating phases are gttasi-Bragg peaks, and 
therefore the corresponding density correlation functions 



display real-space power-law decay, similar to Eq. 1.12 
Although, in principle, the critical values pf, pf m for the 
appearance of each of these floating phases are most likely 
distinct, for simplicity of the presentation we have taken 
Pc = Pc" 1 = Pc- If m reality these critical values are 
sufficiently distinct, and p® < p£ m , then we expect an 
intermediate range of p values, pf < p < p£ m , for which 
no FSm appears. 
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are in the same universality class as the well-known ther- 
mal roughening transition, (35| the dual of the Kosterlitz- 
Thouless transition, with an identical phenomenology. 
Similar to the XY-melting LFS-LSm transition discussed 
above, the melting of the FS into the FSm proceeds via 
unbinding of dislocation pairs with ^-directed Burgers 
vectors. However, because of the presence of massless 
spectator phonon modes in y-direction (transverse to the 
troughs of the periodic potential) , the melting of the FS 
into the FSm is in a new universality class. 

The direct transition from the LFS to the FSm phase 
is most likely first order. Here the order of the u^-modes 
changes from quasi-long-range to short-range (via un- 
binding of type I dislocations, i.e., those with Burgers 
vector parallel to the troughs of the periodic potential.) 
and for the Uy-modes from long-range to quasi-long-range 
(via depinning from the laser potential, i.e., a roughen- 
ing transition) . If both order parameters become critical 
at the same point in the phase diagram, which will be 
the case at the multicritical points where the FS, LFS, 
and FSm phases meet, we have a phase transition corre- 
sponds a simultaneous transition of the KT type and its 
dual analog. 

The remainder of this paper is organized as follows: 
in Section ||, we introduce and motivate our model for 
2d solids subjected to a Id periodic potential and dis- 
cuss the details specific to the experiments on 2d colloids 
in the laser potential. [^2| In Sec. Ill wc give a detailed 
analysis of all the phases which are allowed by symmetry. 
In particular, the static structure factors and correlation 
functions are discussed. The mechanisms, dislocation un- 
binding and soliton proliferation, driving the phase tran- 
sitions are investigated in Sec. IV. In Sec.|v| we derive the 
universal features of the melting phase boundary, demon- 
strating that for sufficiently short-range interactions it 
generically exhibits a reentrant melting observed in the 
experiments of Wei etal. . |2|] Some aspects of the re- 
sponse of the translational and bond-orientational order 
parameter to a small external Id periodic potential are 
analyzed in Sec. VI using a renormlization group crossover 
analysis. In Sec. VII we elaborate on some implications of 



u 



our results to experiments and for computer simulations. 



II. BASIC INGREDIENTS 



K 



FIG. 5. Schematic phase diagram for a primary commen- 
surate orientation with commensurability parameter p > p c 
(The case p = 4 is shown here). As in Fig. |^ the thick line 
indicates a first order transition. Insets: Schematic structure 
factors. As in Fig.^, the x's indicate delta-function Bragg 
peaks and the shaded circles algebraic peaks. 

Wc find that phase transitions between the correspond- 
ing locked and floating phases (LFS-FS and LSm-FSm) 



A. "Microscopic" model 



In the absence of external perturbations, we expect 
that, at sufficiently low temperatures the 2d colloidal 
system freezes into a hexagonal 2d crystal illustrated in 
Fig-! 
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H n = - 



FIG. 6. Perturbation-free ideal hexagonal colloidal crystal, 
characterized by fundamental lattice vectors e^. 

Its lattice sites r„ = niei + 7i2e2, with 774.2 € 2 are 
spanned by a set of lattice vectors 



ei = ae x , 
a 



e2 = -^(e x + V3e y ) , 
e3 = - V3e y ) . 



(2.1a) 
(2.1b) 

(2.1c) 



or equivalently, in Fourier space, the lattice is charac- 
terized by a set of three fundamental reciprocal lattice 
vectors 



Gi = (2n/a')e y , 

G 2 = {7r/a')(V3e x - e y ) , 

G 3 = {n/a'){V3e x + e v ) , 



(2.2a) 
(2.2b) 
(2.2c) 



with a' — and a the mean colloidal spacing re- 

lated to the particle density p by p — 2/\/3a 2 . 




FIG. 7. A set of three fundamental reciprocal lattice vec- 
tors Gi, which completely characterize a perfect hexagonal 
lattice. 

At sufficiently long scales |3(| and to quadratic order 
in the elastic strain 



1 



-(diUj 



djUi) 



(2.3) 



associated with the colloidal displacement field u(x, y) 
relative to the equilibrium position in the unconstrained 
solid, the elastic energy of a 2d hexagonal crystal is well 
described by the continuum isotropic elastic Hamiltonian 



d 2 r 



(2pu 



(2.4) 



The Lame coefficients \x and A, with // the usual shear 
modulus are the only two elastic constants necessary to 
completely characterized the elastic energy associated 
with small deformations of an unperturbed 2d hexago- 
nal solid. 

An applied Id periodic potential, which in experiments 
Jl6| , p2] | with dielectric colloidal spheres is conveniently 
created by two interfering laser beams, is easily incorpo- 
rated as an additional energetic contribution 

ffK = -E/R-^^cos[K-(r„ + u(r„))] , (2.5) 



where we have focussed on the energetically most impor- 
tant lowest harmonic K of such a laser-induced potential. 
The coupling Uk measures the strength of the Kth har- 
monic of the laser potential, proportional to the input 
laser intensity 

For the purpose of the discussion in this section we 
have chosen an "internal" reference frame (e x ,e y ) where 
the orientation of the hexagonal lattice is kept fixed. 
Later, beginning with Sec. Ill we will switch to a "lab- 
oratory frame" (x, y) where the orientation of the laser 
potential is fixed with K || y, i.e. the troughs are running 
parallel to the x-axis. 



B. Commensurability and reciprocal lattice 

For a general wave vector K, the periodic (laser) po- 
tential — characterized by a plane wave, e zK r — will 
not be commensurate with the hexagonal lattice. Only 
for a particular orientation and magnitude of K will the 
spacing between the potential minima match with the pe- 
riodicity of the hexagonal lattice. It is this special set of 
commensurate periodic potentials that is the focus of our 
work here. The characteristic set of commensurate wave 
vectors is easy to find since the reciprocal lattice is defined 
to be the set of all wave vectors G that yield plane waves 
with the periodicity of a given Bravais lattice. Hence, 
commensurability is equivalent to the condition that K 
coincides with one of the reciprocal lattice vectors G. 

In other words, the planes defined by the minima of 
the external potential (cos(K • r)) are a superset of the 
family of lattice planes (Bragg planes) defined by the 
shortest reciprocal lattice vector, G,^ = miGi + TO2G2 
with Miller indices mi and m2 , parallel to the wave vec- 
tor of the external potential K, 



K = pG.,% = piGi + p 2 G 2 



(2.6) 



Note that here we focus on situations where the colloidal 
particles are allowed to sit in the minima of the external 
potential only. More generally, one could also consider 
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situations where the commensurability parameter is less 
than 1 with p a rational fraction. |3(i[| 

With d = 2n/ | K | being the periodicity of the poten- 
tial and a'^ = 2tt/ \ | defining the distance between 
the lattice planes, the commensurability ratio p is given 
by 



K 



d | 
V3a/2 



G" rn. 



m 1 m 2 ) 



-1/2 



(2.7a) 
(2.7b) 



This allows us to characterize the laser potential by a set 
of Miller indices m = (mi,mj) and the commensurabil- 
ity ratio p, i.e., in summary by a commensurability vector 
p = pfh = {pi,P2)- Equivalently, the orientation of the 
Bragg planes can also be characterized by the shortest 
direct lattice vector pointing parallel to the troughs of 
the external potential, 



(2.8) 




FIG. 9. 2d hexagonal colloidal crystal in the presence of a 
commensurate Id periodic potential with period d, commen- 
surability vector p— (1, —1), and potential maxima indicated 
by solid lines. Dashed lines denote the Bragg rows picked out 
by the laser potential minima. 



with the condition R^j • G,^ = 0, i.e., (ni,7i2) = 
(mi, —m<i) a set of integers (direct lattice Miller indices), 
with no common factor complementary to the Miller in- 
dices. 

In Figjj] we displayed an example for the simplest 
set of relative orientations between the periodic poten- 
tial and the colloidal crystal. In our notation it corre- 
sponds to an orientation (mi, 7712) = (1,0) (or equiva- 
lently (ni,ri2) = (1)0)) and a commensurability ratio 
p = 3, i.e., K = 3G(i.o) = 3Gi and a Bragg row spacing 
a,'- = a' = a^/3/2. Because in such m = (1,0) orienta- 
tions, it is primary Bragg rows p8| that run parallel to 
the periodic potential troughs, we call these relative ori- 
entations "primary" . Aside from the simplicity of these 
configurations, our interest in them is driven by experi- 
ments in Refs. |l6| , p2| , where a primary p — 1 orientation 
was studied. 



In addition to these primary p = p(l, 0) configurations, 
we will also make detailed predictions for the next sim- 
plest p = p(l,—l) set of relative lattice-laser potential 
configurations, illustrated for p = 1 in Fig.^. We call 
these orientations "dual-primary" , because they corre- 
spond to Bragg rows running perpendicular to a funda- 
mental real space lattice vector with K = p(Gi — G2) = 
— e^An/a 2 , rather than to one of the three fundamental 
reciprocal lattice vectors. In terms of the direct lattice 
these dual-primary orientations correspond to (n\, n^) = 
(1,1) and Bragg row spacing a' ifi = a/2. 



m • * 1- 




FIG. 8. 2d hexagonal colloidal crystal in the presence of a 
commensurate Id periodic potential with period d, commen- 
surability vector p — 3(1, 0), and potential maxima indicated 
by solid lines. Dashed lines denote the Bragg rows picked out 
by the laser potential minima. 



FIG. 10. 2d hexagonal colloidal crystal in the presence of a 
commensurate Id periodic potential with period d, commen- 
surability vector p— (2,-1), and potential maxima indicated 
by solid lines. Dashed lines denote the Bragg rows picked out 
by the laser potential minima. 
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Using the defin ition of commensurate configurations, 
Eq.Ejt in Eq.2.5, we find that Hk_ reduces to 



AT 



(2.9a) 



-UkoT 2 J d 2 rcos[K-u(r)] , (2.9b) 



where in going from Eq. 2.9a to Eq. 2.91; we went over 
to a continuum description, an innocuous approximation 
for the smooth, |u(r n+ i) — u(r„)| -C a distortions, that 
we study here. 

An equivalent "soft-spin" continuum description of 
the above interaction is in terms of the elementary 
translational order parameters yOG;(r) = p°,.e lGi '^ r+u ' r ^ 
with i — 1,2. The laser-induced periodic potential 
hn(r) = Re[h^e lKr ], acts like an ordering field on 
the p = ±(p!,p 2 )th harmonics of the fundamental or- 
der parameters pc (r), allowing a linear coupling to 

H K = -aJ d 2 r [h* K (r)p G (r) + c.c] (2.10a) 



ap G h° K I <Tr 



e i(G-K).r e iG.u + c ^ 



(2.10b) 



which is finite at long scales only if the condition Eq j2.G| is 
satisfied, in which case it reduces to the expression given 
in Eq. 2.9b, with the identification Uk/o, 2 — 2ap G h^. 
Hence the periodic potential explicitly breaks transla- 
tional symmetry and therefore induces a finite transla- 
tional order parameters P±kM throughout the phase di- 
agram. 



C. Broken rotational symmetry and anisotropic 
elasticity 

The imposed Id periodic potential also explicitly 
breaks, throughout our phase diagram, the 2d rotational 
symmetry down to Z2 (Ising) symmetry, corresponding 
to rotations by 7r. We can see this more explicitly by not- 
ing that the laser potential /ikOO = Re [/!^e' K ' r ] gener- 
ates a set of even-rank tensor fields, 

= d n h K (r)d l2 h K (r)...d l2n h K (v) , (2.11) 

where in above the overline denotes a spatial average. 
The lowest order, rank 2, tensor field is given by 



(2.12a) 
(2.12b) 



It is clear from their definition, that these laser-generated 
2n-rank tensor fields have strengths proportional to 
(Uk) 2 ™ oc (/in) 2 ". They act as external ordering fields, 



which explicitly break rotational invariance (modulo ro- 
tations by 7r) of our system and therefore induce through- 
out our phase diagram finite 2n-adic orientational order 
parameters. These can be characterized by rank 2n sym- 
metric traceless tensors, which are real irreducible rep- 
resentations of the rotation group and are the "angular 
harmonics" of the lowest order, rank 2, nematic order 
parameter 



Qif = Sihifhj - ^6 i: 



(2.13) 



The unit vector n, defines the principle axis of the ne- 
matic order, and, given Eq.2.12b, points parallel or per- 



pendicular (depending on the sign of the coupling be- 

(2) (2) 

tween h\j' and Q\j ) to the periodic potential wavevec- 
tor K. In two dimensions these 2n-rank tensor repre- 
sentations are well-known to be equivalent to the one- 
dimensional complex irreducible representations 



lp2n = e 



i2nU 



(2.14) 



Since in the presence of these laser-induced ordering fields 
all ip2n orientational order parameters are finite through- 
out our phase diagram, no sharp continuous orientational 
ordering phase transitions are possible in our system. 
This is in contrast to the melting of the unperturbed 
lattice, where a thermodynamically sharp orientational 
phase transition is allowed between the isotropic and the 
anisotropic (e.g, hexatic, in a hexagonal lattice) liquids. 
H Therefore, throughout this paper we confine our at- 
tention only to phases and phase transitions that sponta- 
neously break the translational symmetry of the explicitly 
orientationally ordered, modulated colloidal liquid phase. 

The existence of these orientational ordering fields hm 
has important consequences to the form of the colloidal 
crystal elastic energy. To deduce the form of the appro- 
priate elastic Hamiltonian it is instructive to first con- 
sider a 2d hexagonal lattice in the absence such explicit 
symmetry breaking fields. Such state is characterized by 
a finite value of the hexagonal orientational order pa- 
rameter tpQ , J5| with the full 2d rotational symmetry bro- 
ken down to the symmetry of discrete rotations by 2tt/6. 
Nevertheless to a quadratic order in the strain tensor Ujj , 
the energy is invariant under a full 2d rotation group. 

In the absence of a periodic potential the hexagonal 
orientational order can be further spontaneously broken 
down to a lower symmetry. A physically important ex- 
ample is a uniaxially distorted hexagonal 2d crystal of 
anisotropic, orientationally ordered molecules, as, for in- 
stance found in a nematic liquid crystal. Such a system 



exhibits a spontaneous nematic order parameter Q 



(2) 



which modifies the isotropic elasticity Hq, Eq. 2.4. To 
a quadratic order in the strain mj three additional ener- 
getic contributions 



SHn = 



d 2 r 



(2) 



>(2) 



(2.15) 
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are allowed. Because the nematic order is induced spon- 
taneously, simultaneous rotations of the lattice degrees 
of freedom and of the nematic axis (encoded in Q|| ), 
relative to an arbitrary frame fixed in the lab, is clearly 
a symmetry of such uniaxially distorted lattice. It is not 
difficult to show that this rotational freedom allows us to 
eliminate a% coupling linear in it^- 



H ai = aiS J 



d 2 r 



1 



sin w - 



u xy 2sin9 cost 



1 



(2.16) 



by a judicious choice of the rotation angle 8 and a uni- 
axial area-preserving distortion 

Ui — > Ui + fa . (2-17) 

It is important to note that this is only possible because 
the three independent degrees of freedom, 8, fa and <f> v , 
at our disposal, are sufficient to cancel the three inde- 



l yy> 



and 



■'■V 



in H a ,, Eq.2.16 



contributions, Eq.2.15, to the 



pendent linear terms u xa 
Adding H a2 and H a3 
Hamiltonian of an undistorted hexagonal lattice, Eq. 2.4 



we find a general elastic Hamiltonian for a spontaneously 
uniaxially distorted he xag onal lattice (in the absence of 
an external potential) |B7j is given by 



2^L 



uy 



A 



xy ^xx ^yy 



(2.18) 



where we have chosen a coordinate system in which the 
x and y axes coincide with the n and z x n principle 

(2) 

axes of the orientational nematic order parameter Q\- . 
The two additional elastic constants, a total of four in 
Hq , are consistent with two new couplings cti and al- 
lowed by the finite orientational nematic order parameter 

(2) 

Qij ■ The four independent elastic constants also coin- 
cide with the expectation, that with a symmetry between 
x and y broken, the elastic energies associated with the 
strain tensor components 



*yy ~ "y^y > 



(2.19a) 
(2.19b) 

(2.19c) 



are clearly independent. Although rotations relative 
to the orientational uniaxial order is no longer a sym- 
metry of Hq, because only the symmetric strain tensor 
Uij enters the elastic energy, H§ is still invariant under 
"atomic" displacements 



u = vz x r 



which correspond to global rigid rotations of the 2d solid, 
by an infinitesimal angle 8 about the z-axis. This lat- 
ter symmetry is clearly present in an anisotropic lattice 
without an external pinning potential. 

In contrast, however, in our system, the Id periodic 
potential has a fixed orientation in the laboratory frame. 
Hence, in addition to the uniaxial lattice anisotropy, such 
a potential also explicitly breaks symmetry of rotations 
relative to the lab (and therefore to the periodic poten- 
tial) frame. It therefore picks out a special coordinate 
system relative to which the angle 8 is measured. 

As discussed above such external potential acts as an 
external 2n-rank tensor fields and explicitly breaks the 
corresponding orientational symmetry. The appropriate 
elastic energy can be deduced by focusing on the lowest 

(2) 

order nematic ordering field h\j . It allows the following 
additional energetic contributions 



Hh 2 



d 2 r 



(2.21) 



that explicitly break symmetry of rotations relative to 
the frame picked out by the periodic potential. 

For the purposes of classification of the rel ative orien- 
tations discussed in previous subsection, Sec. [IB , it was 
more convenient to keep the lattice fixed and to rotate the 
periodic potential into a particular orientation, uniquely 
labeled by an integer 2d Miller index vector p — {p\,P2)- 
However, once an orientation has been selected and clas- 
sified by p, to analyze the continuum elastic model and 
its thermodynamics that follows it is more convenient 
to work in a coordinate system in which, instead, the 
troughs of the Id periodic potential run along the new 
x-axis. For such a choice of a lab coordinate system, 



,( 2 ) - 



(2.22) 



(2) 

Using this expr essions for together with Q.^ , Eq 



inside Eq. 2.21, and combining it with H ai , Eq.2.16, we 



2.13 



find the following symmetry breaking energetic contribu- 
tion, which, in the presence of a Id periodic potential 



must be added to Hq, Eq. 2.18 



H 



ai+h 2 



d 2 r 



a!Su xx (sm 



) + u yy(0ilS COS 2 



ex S 1 
— h) + a 1 Su xy sin28 - -hScos28 , (2.23) 



where h = ^ 



K 2 1 h | 2 and angle 9 measures the deviation 
of the nematic axis n away from K set by the orienta- 
tion of the periodic potential. While it is still possible 
to eliminate the terms linear in u xx and u vy by a lattice 



distortion Eq.2.17, in the presence of the external poten- 
tial it is no longer possible to shift away the u xy term. 
Selecting fa so as to cancel u XXl u yy and combining the 



(2.20) resulting H ai+ h 2 with Hq we find 
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H a = / d 2 r 



X X yU xx u yy + au xy sin 20 - - cos 20 



(2.24) 



where we defined rotational symmetry breaking couplings 
a = a\S and 7 = hS. It is clear from the above a and 
7 terms in H a that, in the absence of strain, u xy = 0, 
the energy is minimized by = 0, corresponding to the 
nematic axis alignment with K, imposed by the periodic 
potential. In the presence of fluctuations will be small 
but finite. Expanding H a , above, in these small fluctu- 
ations, we obtain a final form of the elastic Hamiltonian 
characterizing our system 



Ho: 



d 2 r 



1- , 1 



%l lu X y + 2^ xxU xx + 2^ yyU yy 



KyU xx u yy + 2au xy + 2j0* 



(2.25) 




FIG. 11. A u xy 7^ shear deformation (shown for simplic- 
ity for a square lattice) with principle axes along the (1, 1) 
and (1,-1) directions in the xy-plane. In the presence of a 
trough potential (dashed lines) parallel to the x-direction, the 
particle array, with axes x' and y' , can lower its interaction 
energy with the periodic potential by rotating in a clockwise 
direction to bring the particles into better alignment with the 
minima in the trough potential. 



Now to complete our derivation we must relate the an- 
gle 0, characterizing the orientation of the nematic order 
to the elastic Ui degrees of freedom. We expect that the 
orientations of the nematic and hexatic order parameters, 
present in our uniaxially distorted hexagonal lattice, are 
locked together. Since in the crystalline phase, fluctua- 
tions in this bond orientational order are in turn locked 
to the local rotation angle induce d by the phonon dis- 
placements, in the Hamiltonian Eq. 2.25| we can make the 
well-known identification 



Eq.pT25 



together with 
Eq.pbl, de- 



The elastic Hamiltonian H e \ 
the commensurate pinning potential i?K 
fines our model 2d colloidal system in the presence of a 
commensurate periodic laser potential. Our aim, in the 
remainder of the paper is to analyze the symmetry al- 
lowed phases and the nature of the transitions between 
them embodied in this model. 



III. SYMMETRY- ALLOWED PHASES 



(8 X Uy - 8yU X ) 



(2.26) 



thereby completing our derivation. We find that the re- 
sulting elastic Hamiltonian, which characterizes a hexag- 
onal lattice in the presence of a Id periodic potential, in- 
volves 6 elastic constants. While a similar form was sug- 
gested, based on symmetry, by Ostlund and Halperin BTf] 
in their analysis of melting of distorted hexagonal cr ystal 
films, the a term appearing in our Hamiltonian, Eq. 2.25 
was missed in their expression. As illustrated in Fig.|ll|, 
physically, the a term is present because, with troughs 
running along the x-direction, an xy-strain will bring par- 
ticles in Bragg planes lying in the troughs out of align- 
ment with the minima of the periodic potential. This 
generates a torque which attempts to rotate the lattice 
and improve the alignment. 



The starting point of our analysis is the model Hamil- 
toni an H = H e \ + Hj£, obtained from combining Eqs.2.9b 
and [2.2li . Here we have chosen (without loss of general- 
ity) K to lie along the y-axis, i.e., the periodic potential 
troughs running parallel to the convention that 

we will stick to throughout the remainder of the paper. 
This Hamiltonian admits a rich variety of thermodynam- 
ically distinct phases. As discussed in the Introduction, 
the phase diagrams depend on the commensurability ra- 
tio p, or more specifically, in which of the three regimes 
p = 1, 1 < p < p c , or p > p c , p actually falls. The com- 
plexity of the phase diagram is highest for p > p c , and so 
in order to discuss all the phases possible in our system, 
we focus on these high p commensurability ratios. 

It is convenient to enumerate the five allowed phases 
starting with the most ordered, which naturally occurs 
at the lowest temperatures, and proceeding toward the 
higher temperature disordered phases, by invoking two 
types of disordering mechanisms, dislocation unbinding 
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and soliton proliferation. A detailed investigation of 
these mechanisms is deferred to the subsequent section, 
where we discuss the nature of these transitions and their 
hierarchy as a function of temperature and periodic po- 
tential strength (laser intensity). Here, we focus on the 
phases themselves, rather than on their location in the 
resulting phase diagram. As discussed in detail in the 
previous section and in the Introduction, the imposed 
periodic potential explicitly breaks rotational symmetry, 
and therefore all five phases exhibit true long-range ori- 
entational order. This external potential also explicitly 
breaks continuous translational symmetry along y (with 
potential troughs taken to run along x) down to a discrete 
symmetry of translations by the period d of the potential. 
Hence all phases will trivially exhibit long-range order in 
the corresponding translational order parameter, leading 
to true delta-function Bragg peaks at the multiples of 
the reciprocal lattice vector (27r/rf)y in their structure 
functions. 



A. Solid Phases 



As in the absence of a periodic potential, the most 
ordered phase of isotropic particles confined to 2d is a 
solid. The striking effect of turning on an external Id 
periodic potential is that it can lead to two thermody- 
namically distinct uniaxially distorted hexagonal crystal 
phases, which we term a "locked floating solid" (LFS) 
and a "floating solid" (FS). Being crystalline, both of 
these phases exhibit 2d translational (quasi-long-range) 
order, and are characterized by a finite shear modulus. 
These emerge as a result of breaking the translational 
symmetry T% <g> T x of the "modulated liquid" (ML) (see 
below), corresponding to independent discrete transla- 
tions by dy and continuous translations along x, down 
to 2d discre te transla tions generated by lattice vectors ei 
and B2, Eq. 2. la, 2.1b. Although in the presence of ther- 
mal fluctuations dislocations will be thermally nucleated, 
in the solid phases they are confined to finite size dipoles 
with a zero Burgers "charge". These, therefore, can be 
safely integrated out of the partition function with only 
weak finite renormalization of the elastic constants. Con- 
sequently, a purely elastic description in terms of phonon 
modes u x and u y is appropriate in both phases. 

The LFS and FS phases differ in the importance of the 
periodic pinning potential. In the FS, expected to occur 
at temperatures higher than the LFS, thermal fluctua- 
tions in the positions of the colloidal particles are suffi- 
ciently large such that at long length scales they average 
away most J38| of the long scale effects of the periodic 
potential. In contrast, in the LFS the periodic poten- 
tial strongly pins the colloidal crystal transversely to its 
troughs. 



1. Floating Solid (FS) 

The floating solid can be rigorously differentiated from 
its locked counterpart as a 2d colloidal crystal phase in 
which the periodic potential is irrelevant in the renor- 
malization group sense. This implies that at long scales, 
many, but not all (see below) of the thermodynamic prop- 
erties of the FS ar e well described by the elastic Hamilto- 
nian H c \, Eq. 2.25 , with two coupled "massless" u x and u y 
degrees of freedom, and ignoring H^. Therefore, in many 
ways the FS is qualitatively quite similar to a 2d solid 
without the periodic pinning potential. In particular, this 
similarity extends to the lattice displacement correlation 
functions which are logarithmic in x and y. However, 
these similarities do not extend to all correlation func- 
tions, and the periodic pinning potential has important 
qualitative consequences for the FS phase that distin- 
guish it from an ordinary 2d solid. The "irrelevance" 
of the periodic potential means only that a perturbative 
expansion in Uk, for a sufficiently small value is conver- 
gent. Consequently, average quantities, that are finite 
at Uk — can be well approximated by their Uk = 0- 
values, i.e., working with H « H e \, as is usually done. 
However, quantities that vanish (or diverge) to this ze- 
roth order, must be evaluated to the next lowest order in 
Uk to obtain a nontrivial (finite) result. 

To illustrate this point, recall that the periodic poten- 
tial explicitly breaks rotational and translational symme- 
try, despite its irrelevance in the FS phase. While the for- 
mer leads to uniaxial anisotropy in the hexagonal lattice, 
the latter is responsible for the true long-range order in 
the translational order parameter p G (r), with G integer 
multiples of the wave vector K characterizing the peri- 
odicity of the external potential. In the presence of the 
periodic potential, even the most disordered modulated 
liquid phase (see below) displays true long-range trans- 
lational and orientational order. Clearly then, a more 
ordered FS will also break these symmetries. 

As a concrete example of how the periodic potential 
affects the FS phase, consider the real-space 2-point cor- 
relation function of the translational order parameter 



Pair) 



= e iG-u(r) 



defined by 



C G (r) = (p G (r)p G (0)) , 

= / e »G.(u(r)-u(0))j ( 

^c£(r) + {p G )(p G ) 



(3.1) 

(3.2a) 
(3.2b) 
(3.2c) 



where in Eq.3.2c, C^T^r) is the connected part of C G (r). 



The distinguishing feature of the FS phase is the irrele- 
vance of the periodic potential i?K- Hence in the limit 
of a weak laser potential, i.e., small Uk, we can compute 
C G (r) in a controlled, convergent perturbative expansion 

in Uk- The connected part Cq (r) is nontrivial even to 
zeroth order in Uk , and a standard calculation gives, || 
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Mia 



where 



|G| 2 k B T 3/i + A 

4"7T [I 2p, + A 



(3.3) 



(3.4) 



For simplicity, we have used the isotropic elastic Hamil- 
tonian, Eq.2.4, in place of the cor rect s ix elastic constant 
anisotropic Hamiltonian H c \, Eq. [2.25 , which leads to a 
qualitatively similar, but anisotropic power-law decay of 
spatial correlations. We use long wavelength elastic con- 
stants finitely renormalized by thermally excited bound 
dislocation dipoles. 



We can compute the persistent part of Cg (r) , by calcu- 
lating 



(Pa) = 



iG-u(°)\ 



(3.5) 



in a perturbation theory in Uk , which, because of the ir- 
relevance of the periodic potential is convergent in the 
FS phase. For Uk — the translational order pa- 
rameter vanishes like (pg) = {L/a)^ r,a ^ 2 with system 
size L — * oo. Upon expanding the Boltzmann weight 
e -(H +H K )/k B T j n a p 0wer series in Uk, we find to lead- 
ing order in Uk and for L — > oo 



(p- 



G/ 



OO ^ 

n=l 



U K 



2a 2 k B T 



ft n lb '11 

(*G,nK + *G,-nK) R / d 2 r 3 gflrjl/ a)^ JJ(|r, - Tjl/a)^* 

A 1 w A 1 i ^ A 



3=1' 



i<3 



(3.6) 



Here we have used the fact that for L — > oo 



(exp 



u(r Q ) 



exp 



a</3 



q/3 



G {c \v a 



(3.7) 



for J2 a la — and zero otherwise. We have also intro- 
duced a phonon connected correlation function G^ (r) 



G( c >(r)^(|u(0)-u(r)| 2 ) o 



(3.8) 



Averages with elastic Hamiltonian are designated by 
(...)o- Upon again approximating H e \ by its isotropic 
form Ho, Eq.2.4, a straightforward calculation in the 
limit L/a^S> 1, r/a» 1, gives 



G (c) (r)Rj g| log(r/c 



(3.9) 



Since Uk is irrelevant in the floating solid phase the in- 
tegrals in Eq.3.6 are IR convergent (i.e. for L — > oo). 
The power-laws appearing in the integrand are implicitly 
understood to be cutoff below the lattice constant a scale 
by the obvious behavior (see Eq.|3.8D of the phonon cor- 
relation function lim r _> Q (r) = 0. Upon performing 
the spatial integrals, which are dominated by the behav- 
ior of the connected phonon correlation function at small 
distances (UV, lattice cutoff a), we obtain up to non- 
universal factors of order 1 



(p, 



1 



G 



U K 



2a 2 k B T 



(*G,nK + *G,-nK) • (3-10) 



Hence, as argued above on physical grounds, despite 
of the irrelevance of the periodic potential, the FS dis- 
plays true long-range order in the translational order pa- 
rameter pa, with G satisfying G = ±nK, with C(j(r) 
approaching its asymptotic value as a power-law in r. 



Other translational order parameters, with G not satisfy- 
ing the above condition have pure power-law correlation 
functions, decaying to zero at long separations. In par- 
ticular, these include the fundamental translational order 
parameters pGi , which display quasi-long-range order in 
the FS phase. 

Having calculated the translational correlation func- 
tion Cg(i"), the structure function 



S(q) 



-Zq-r„ 



Cq( r r, 



(3.11) 



can now be easily obtained. Using Eqs. 3.2a, 3.3, and 3.10 



and taking advantage of the Poisson summation formula 
to perform the sum over the lattice sites r m , we find 



S(q) » £ 



1 



G 



q- G| 2 "^G 



71— — OO 



with 



i / u 



K 



(n!) 2 \2a 2 k B T 



2 ii 



(3.12) 



(3.13) 



and prime on the summation in Eq.3.12 indicating that 
the n = term is excluded. 



Equation 3.12 predicts true Bragg peaks (with power- 
law corrections) at multiples of the periodic potential 
wave vector K and pure power-law (quasi-) Bragg peaks 
at all other reciprocal lattice vectors G, even for those 
with G || K. Note that in a real physical system, the pe- 
riodic potential will not in general be a single harmonic as 
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assumed in our model, Eq.2.5. Hence, we expect that the 
Bragg peak amplitude A n observ ed in experiments will 
be a sum of terms like those in Eq. 3.13 and the square of 



the amplitude of the nth Fourier harmonic, £/ n K of the 
applied periodic potential. This of course will only mod- 
ify the prefactors in the di fferen t contributions to 5(q), 
predicted for the FS in Eq.3.12. We schematically illus- 
trate 5(q) for a floating solid in Fig.|l2| for the commen- 
surability vectors p = (5,0) andp = (2, —2), respectively, 
with the y-axis chosen to point along K. 



(a) 



-e—)t 



(b) 



« 

o q x 



FIG. 12. Schematic structure function for the FS phase 
with the commensurability vector (a) p = (5, 0) and (b) 
p = (2, —2), illustrating a comb ination of the quasi- and true 



Bragg peaks, given by Eq.3.11, Crosses indicate true Bragg 



peaks and open circles quasi-Bragg peaks. 



pinning energy Ujc is much larger than the elastic en- 
ergy /j,a 2 for the shortest (and therefore all) wavelength 
phonon mode, our system is in the strong pinning regime. 
For a commensurate periodic potential, in this regime, 
fluctuations in the lattice positions perpendicular to the 



troughs are small and the periodic potential Hj^, Eq. 2.9 



can be safely expanded in powers of the corresponding 
phonon degree of freedom, K • u, leading to 



Hi 



K 



const. 



2 UkC1 
1 

2 l 



d 2 r(K- u(r) 



const. + ^U K a 2 K 2 j d 2 ru 2 y {Y) 



(3.14a) 
(3.14b) 



In contrast, a weak pinning regime, Uk <C /ia 2 , consists 
of two sets of elastic modes, those with k < k c and those 
with k > k c , where k c = K/b* is a crossover wavevector 
for which the elastic energy density n(k c a) 2 just balances 
the pinning energy density Uk^^K 2 at the same length 
scale. Since the pinning energy is subdominate to the 
elastic energy for modes with k > k c , we can simply in- 
tegrate out these weakly pinned modes perturbatively in 
Uk- This results in an effective strength of the pinning 
potential given by 



U K (K) = U K b, 



-W2 



(3.15) 



After equating this to the corresponding elastic energy 
/i(a/b*) 2 we find 



6* 



/ / 1! 

Uk~ 



2\ 2/(4-»7 K ) 



which, when inserted inside Eq.3.15 leads to 

W(4-T7k) 



EMM = Uk 



\fJM 2 



^4/(4-^)^^2^/(4-^) 



(3.16) 

(3.17a) 
(3.17b) 



The set of on-q^-axis quasi-Bragg peaks (indicated by 
open circles) interleaving the true Bragg peaks (indicated 
by "x" 's) is the notable feature that distinguishes the FS 
from its locked counterpart LFS, in which all on-qj,-axis 
peaks are true Bragg peaks. 



2. Locked Floating Solid (LFS) 

At sufficiently low temperatures, the periodic poten- 
tial will always be a relevant perturbation, pinning the 
2d solid in the direction perpendicular to its troughs. Be- 
cause of the Id nature of the pinning potential, the 2d 
crystal will remain unpinned along the direction of the 
potential minima and will be able to adjust freely in that 
direction. To reflect this dual character, we therefore call 
this phase the "locked floating solid". 

At high laser intensity, such that the bare value of the 



Since the u y fluctuations in the remaining strongly 
pinned elastic modes are small, the effective pinning po- 
tential i?K can once again be safely expanded in powers 



of u y . Doing this we obtain a result identical to Eq.3.14b 



but with Uk replaced by Uk(}>*) given in Eq.3.17. 

Hence, in both the strongly and weakly pinned regimes, 
unlike the FS phase, the LFS is characterized at long 
wavelengths by one acoustic (u x ) and one optical (u y ) 
phonon mode, with an effective Hamiltonian 



H = H cl + ^J d 2 ru 2 y (v) 



(3.18) 



Here, we have introduced a correlation length £ which, 
given Eqs.3.14t, 3.17, reads 



C 2 (Uk) 



^K 2 , 

fia z 



^ ^ » 1, 



K 2 , for ^ < 1. 



(3.19) 
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At lengt h scales longer than the crossover scale set by 
£, Eg. 3. 19 , we can safely ignore the spatial derivative of 
U y terms, and the LFS is well described by an effective 
Hamiltonian 



#LFS = ^ 



where 



d 2 r 



P_ 2 

£2 u y 



B yx = {p + 7 - a) , 

Bxx = ^XX • 



(3.20) 



(3.21) 
(3.22) 



We can now compute the translational order parame- 
ter correlation function and the structure function that 
characterize the LFS phase. Repeating first the calcu- 
lation for the persistent part determined by (pc), we 
immediately find, that, as in all the phases in the pres- 
ence of the periodic potential, (pc) ^ for G = ±nK. 
However, the distinguishing feature of the LFS is that 
this average is finite for all G parallel to K, by virtue 
of the finite pinning length £, Eq.3.19. This result can 
be immediately seen by noting that for G||K, the loga- 
rithmically divergent (with L) (u x )o correlation function 
does not appear in (pc)o, where the subscript "0" again 
represents an average with the elastic Hamiltonian H e i 
only. Instead we have 



G 



- P -hG 2 (ui) 



Va/2 



(3.23a) 
(3.23b) 



which only involves the "massive" u y degree of freedom, 
whose logarithmic correlations are cutoff at L c — £ and 
therefore is finite even in the thermodynamic limit. 

We can also obtain the above result via a straightfor- 
ward matching calculation. The difficulty of computing 
translational correlation functions in the weakly pinned 
regime of the LFS phase is that for long length scales 
(> £), despite of the weakness of the pinning potential, a 
direct perturbative expansion in Uk is divergent because 
of its relevance (in the renormalization group sense) in- 
side the LFS phase. The power of the renormalization 
group is that it allows us to relate this difficult weakly 
pinned, small Uk regime to the strongly pinned regime, 
where Uk has grown to the magnitude of the elastic en- 
erg y pa 2 , and can therefore be treated as a "mass" , as in 
Eq. |3.141: . We can apply this matching procedure to the 
computation of (pg{Uk)) , by using a relation between 
the weakly and strongly pinned regimes, namely 

(pg(Uk)) = b-^ 2 (p G (U K b 2 -^/ 2 )) (3.24) 

obtained using the scaling dimension of the operator p G 
an d th e RG eigenvalue of Uk, both easily extracted from 
Eq.|3.3|. Choosing the arbitrary rescaling factor b = 6„ 



such that Uk (b) i s in t he strongly pinned regime, where 



UK{b*) = pa , Eq. |3.24| becomes 



(pg(U k )) 



pa z 



(p G (pa 2 )) . (3.25) 



Since the right hand side is in the strong coupling regime, 
it can be ea sily c omputed using the coarse-grained Hamil- 
tonian, Eq. 3.18. Doing this we find 



(pg(Uk)) = 



Uk 
pa z 



e 2 



Vg In(Ko) 



(3.26) 



which in the wea kly pin ned regime is equivalent to the 
result given in Eq. 3.23b . 

Note that the nontrivial nonlinear power-law response 
of the translational order parameter to the periodic laser 
potential, predicted by Eq, 3.26| is only a nonanalytic 
piece of the full response, which includes an analytical 
background. Hence, although at low temperatures, such 
that ?7g/(4 — Vk) < L the full response in the Uk — > 



limit is dominated by the nonanalytical part, Eq.3.26 



at higher temperatures, the ever-present linear piece of 
the analytical part will dominate, and experimentally one 
should instead observe 



(pg(U k )) - U K 



(3.27) 



For our highly anisotropic system, the connected part 
of the correlation function Cg(i"), is given by 

C { £\r) = e -k[GlG xx {T)+2G x G y G*y{r)+GlG yy (r)\ ^ (3 28 ) 

where (r) = ((ui(r)—Ui(0))(uj(r)— Uj(0))), is the con- 
nected phonon correlation function computed with the 
full Hamiltonian. In the weakly pinned regime, for small 
length scales, all phonon correlation functions display the 
usual 2d logarithmic growth, which, in the isotropic ap- 
proximation, i.e., using Hamiltonian Hq, Eq.2.4 leads to 



the po wer-law correlation for (r) that we found in 
Eq.3.3 for t he FS phase. However, for length scales longer 
than £, Eq. 3. IE, while G xx (v) will continue to grow log- 
arithmically, such growth in G yx (r) and G yy (r) will be 
cutoff by the pinning length £. Consequently, in the LFS 
we find 



G 



(r) - P 
V r 



r lG v 



where 



G„ 



fc R T 



2 71 " y B xx B y 



(3.29) 



(3.30) 



G .' 

which reduces to rin — tt* , B when the effect of 

the periodic potential on elasticity and renormalizations 
due to dislocation pairs on the effective elastic coefficients 
are neglected. A discrete Fourier transform of this corre- 
lation function gives the corresponding structure function 



k B T 
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S(q) a £ 



a 



B G 



|q- G| 2 -^ 



A G <5( 2 )(q-G) 



(3.31) 



where the quasi-Bragg peak amplitude Bq is given by 



Bq cx 



-ria v /2 
K ' 

K ' 



for ^ > 1, 
for % < 1. 



(3.32a) 
(3.32b) 



and the Bragg peak amplitude Aq 




K ) 

2Wg/( 4 -W k ) 



for ^ 
for ^ 



»i, 



(3.33a) 
(3.33b) 



which is finite if and only if G is par allel t o K. As a con- 
sequence of the discussion after Eq. |3.26| , the amplitude 
Aq will also have a background analytic in Uk, which 
in a weak pinning limit scales as U\ (see Eq. |3.27 ), and 
therefore at higher temperatures w ill dom inate over the 
nonanalytical part predicted in Eq. [3.33b| . 

We illustrate schematically S(q) in Fig.[l^ for the com- 
mensurability vectors p = (5, 0) and p = (2, —2), respec- 
tively, with the j/-axis chosen to point along K. 



(a) 



(b) 
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FIG. 13. Schematic structure function for the LFS phase 
with the commensurability vector (a) p = (5, 0) and (b) 
p — (2, —2), illustrating a comb ination of the quasi- and true 
Bragg peaks, given by Eq.3.11. 



These predictions for the structure function of the LFS, 
displaying amplitudes that vanish as nontrivial powers 
(determined by a continuously varying exponent fj Gy ) of 
the periodic potential strength (Eqs. 3.321, [3.33b ) pro- 
vide the first theoretical explanation for observations of 
Clark, et al. pi. 



B. Smectic Phases 

As was first pointed out by Ostlund and Halperin |57| , 
in uniaxial two-dimensional lattices, dislocations with 
Burgers vector along and perpendicular to the uniaxial 
axis will generically have different core energies and will 
therefore proliferate at different temperatures. This will 
consequently allow the possibility of a phase that is inter- 
mediate between a fully ordered crystal and a completely 
disordered liquid. 

In a commensurate orientation, such that Bragg rows 
coincide with the periodic potential troughs, we would 
expect dislocation pairs with Burgers vectors parallel to 
the potential minima, to unbind first. We call the re- 
sulting class of thermodynamically distinct phases, smec- 
tics. Their main common characteristic is that they 
display a finite elastic modulus for shear deformations 
perpendicular to the Burgers vector of unbound disloca- 
tions, but do not resist shear parallel to them, possessing 
only liquid-like correlations between the corresponding 
"atomic" rows. Consequently, such 2d smectics display 
Id periodicity perpendicular to the Burgers vector of un- 
bound dislocations, and, as illustrated in Fig.[l4, can be 
equivalently described as a periodic stack of Id liquids. 

It is important to note, that despite of their name, the 
smectics discussed here are fundamentally distinct from 
the smectic phases found in liquid crystal materials and 
substrate-free smectics discussed in Ref. j37j. The most 
important distinction is that in liquid crystal smectics 
and those without an underlying pinning substrate, the 
orientational symmetry is broken spont aneo usly (uniax- 
ial anisotropy notwithstanding; see Sec.lIC), leading to 
a soft Laplacian-curvature (rather than gradient-tension) 
elasticity, which preserves this underlying symmetry even 
in the smectic phase, where it is nonlinearly realized. 
p9| , |4C| In fact such substrate- free 2d smectics, because 
of the softness of their elasticity, are well known to be 
unstable to thermally-driven unbinding of dislocations, 
and at scales longer than the distance between these 
free dislocations are therefore indistinguishable from a 
nematically-ordered 2d liquid. Q] As was recognized 
by the authors of Ref. |37]], such thermal instability of 
substrate-free 2d lattices precludes the existence of ther- 
modynamically distinct intermediate 2d smectic phase in 
which only one set of Burgers vectors (e.g., along the uni- 
axial direction) unbind. However, in strong contrast to 
those rotationally invariant systems, in 2d lattices stud- 
ied here the periodic (laser) potential explicitly breaks 
rotational symmetry, binding by a linear potential dis- 
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location pairs with Burgers vector having components 
along K. Consequently, such dislocations remain bound 
even when those with Burgers vectors perpendicular to 
K unbind and therefore allow the existence of 2d smec- 
tic phases that are thermodynamically distinct from a 
liquid. 

Deep in such a smectic phase, the u y (r) phonon field, 
which (see FigjlJ) describes local fluctuations in the max- 
ima positions of the Id density wave, is the only remain- 
ing important degree of freedom. The ever-present bound 
dislocation pairs and the density of vacancies and inter- 
stitials are "massive" degrees of freedom. They can be 
easily integrated out, leading only to a finite renormaliza- 
tion of elastic constants for u y deformations, and there- 
fore are unimportant in a static theory. 
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FIG. 14. 2d colloidal smectic phase in the presence of a 
commensurate Id periodic potential with period d, commen- 
surability parameter p — 3, and potential maxima indicated 
by full horizontal lines. Dashed lines denote the maxima in 
the smectic density, which are pinned inside the minima of 
the periodic laser potential. 

In close analogy to the translational order parameter of 
the 2d crystal, the smectic is distinguished from a liquid 
by a finite translational order para mete r pc = e u , but 
with a single (rather than a set, Eq.2.2) reciprocal vector 
G = Gy = (27r/a)y. It is related to the total molecular 
density via a standard relation |39| 



p(r) =Rc[p + e lG V G (r)] 



(3.34) 



where po is the mean density of the smectic. 

Of course, in the presence of a Id periodic potential, 
a smectic is a thermodynamically distinct phase only if 
G = (27r/a)y differs from the wavevector K character- 
izing the external potential and the modulated liquid. 
Commensurate smectics, which we focus on here, are 
equivalently characterized by the ratio of their period a 
to that of the periodic potential d, with commensurabil- 
ity ratio a/d = p £ Z. A p-Smcctic then spontaneously 
breaks the discrete translational symmetry Tj (g> T x of 
the modulated liquid, with its equal occupancy of each 
potential minima down to T% (g)T x , with only every pth 
minima equivalently populated. Clearly then p = 1 is in- 
distinguishable from the fully disordered modulated liq- 
uid. 

Above symmetry considerations uniquely specify the 
Hamiltonian that characterizes the p-Smectic phase 



S m 



d 2 r 



2 (B X y(d X Uy) 2 + Byy(dyUy) 2 



2 

Uko>~ 2 cos(Ku y (r)) 



(3.35) 



which, not surprisingly, is an anisotropic scalar Sine- 
Gordon model in the phonon field u y (r). 

Given the form of the Hamiltonian in Eq. |3.35 , there is 
a close similarity between the properties of the smectic 
and the 2d crystal, studied in the previous section. The 
quantitative differences between these phases are due to 
the distinction between the vector (u = (u x ,u y )) and 
scalar (u y ) nature of elastic degrees of freedom in the 
2d solid and smectic, respectively. More specifically, in 
close analogy to the 2d solid, we find that for a fixed inte- 
ger commensurability ratio p, there exist a low tempera- 
ture "locked" and higher temperature "floating smectic" 
phase. These are distinguished by the importance of the 
periodic pinning potential, which is relevant (in the RG 
sense) in the LSm phase, acting as a "mass" for u y , and 
irrelevant in the FSm phase, where for most static prop- 
erties it can be ignored. 



1. Floating Smectic (FSm) 

In the "floating smectic phase" (FSm) , thermal fluctu- 
ations in the position of the layers are sufficiently large 
that at long length scales they average away many effects 
of the periodic pinning potential. Hence, many of the 
static properties of the FSm phase can be well described 
by the Hamiltonian, Eq. 3.35| , with Uk = 0. However, as 
we discussed in detail in our analysis of the FS phase, 
despite of the RG irrelevance of the periodic potential, 
continuous translational symmetry is still explicitly bro- 
ken by it, which leads to true long-ranged translational 
order in the smectic order parameter pc for Gy at mul- 
tiples of the reciprocal lattice vector Ky, characterizing 
the laser potential. 



Calculations that closely parallel those of Sec. Ill A 1 
for the FS, lead to power-law correlations in the con- 
nected part of the translational two-point correlation 
function 



4 C) (r) 



1 



where 



VFSr 



r \VFBn 



k B TG 2 



2i:^B^yB 



(3.36) 



(3.37) 



mi 



is the exp onen t characterizing the FSm phase, in analogy 
to?7 G , EqQ of the FS. 

The disconnected part of the smectic translational cor- 
relation function is finite only at G = nK (n £ Z). The 
corresponding floating smectic structure function is given 
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by an expression similar to the FS, Eq, 3.12| . The only 
difference is that fj G of the FS is replaced by rypsm of 
the FSm and the summation over G is a sum over inte- 
ger multiples of 2ir/a. Consequently, one expects to see 
sharp peaks only on the q y axis, with power-law peaks at 
G 7^ nK, and true Bragg peaks at G — nK . This FSm 
structure function is schematically displayed in Fig.[l5| 
below. 



2. Locked Smectic (LSm) 

As the temperature is lowered, the periodic potential 
becomes relevant, pinning the smectic layers. The re- 
sulting "locked smectic" (LSm) phase is characterized 
by long-range translational order, and, as illustrated in 
Fig.[l5|, displays true Bragg peaks at all values of the on- 
q a -axis reciprocal lattice vectors G = nlixja. At long 
scales, the effective elastic Hamiltonian that character- 
izes this phase is simply 



u - 1 I* 



(3.38) 



with £ given by Eq. 3. IS 



(a) 



!! 
!! 



(b) 



FIG. 15. (a) Schematic of the structure function for 
the floating smectic phase, characterized by on-f/y-axis 
quasi-Bragg peaks (open circles) and true Bragg peaks 
(crosses), (b) Schematic of the structure function for the 
locked smectic phase, characterized by on-g^-axis true Bragg 
peaks, with small and large crosses indicating spontaneously 
and directly induced translational order. 



3. Modulated Liquid (ML) 

The modulated liquid is the most disordered phase, 
which occurs at highest temperatures and does not spon- 
taneously break any symmetries. It is characterized by 
a vanishing shear modulus, unbound dislocations, ab- 
sence of massless Goldstone modes, and a discrete sym- 
metry of translations along the y-axis by periodic poten- 
tial constant d. The corresponding structure function of 
this explicitly orientationally ordered phase, illustrated 
in Fig.|l6|, is a set of true Bragg peaks at multiples of 
the reciprocal lattice vector K — 2ir/d oi the periodic 
potential. 



FIG. 16. Schematic of the structure function for the mod- 
ulated liquid phase characterized by on-g^-axis true Bragg 
peaks located at n{2ir/d). 

Finite linear translational order parameter susceptibil- 
ity guarantees that the average order parameter is linear 
in the strengt h of t he periodic potential. Therefore, as is 
clear from Eq.1.17 the strength of the Bragg peaks scales 
as a cube of the input laser intensity, proportional to Uk, 
as observed in experiments by Clark, et al. fljjjLpSl , 



IV. PHASE TRANSITIONS 



Phase transitions that take place in our system fall into 
two broad classes: roughening and melting. However, for 
high values of the commensurability ratio p (j> > p c ) these 
classes are mathematically related to each other by the 
duality transformations |4^,|8| , and are both examples of 
the Kosterlitz-Thouless type of transitions, with kinks 
and dislocations unbinding, respectively. For p < p c , 
the roughening transitions are in a different (Ising model 
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and other models with a discrete symmetry) universality 
class. 



eV 



u<(r) 



V(r') , 



(4.3a) 
(4.3b) 



A. Roughening Transitions 

Phase transitions that fall into the roughening transi- 
tion universality class separate a low temperature ordered 
phase, in which a potential Goldstone mode is strongly 
pinned by an external periodic potential, from a quasi- 
long-range ordered phase, in which the periodic potential 
is irrelevant in a renormalization group sense. The locked 
floating solid to floating solid and the locked smectic to 
floating smectic transitions, discussed in Sec. Ill fall into 
this broad universality class, although differ in details 
that we discuss below. 

Despite these small differences the analysis of these 
transitions are quite similar and can be done via 
standard perturbative momentum-shell renormalization 
group (RG) transformation. j4^^|] Since the smooth 
(locked) and rough (floating) phases are distinguished by 
the relevance and irrelevance of the periodic potential, 
respectively, we can find the trans ition temperature by 
analyzing the behavior of H^, Eq.2.5, as a function of 
length scale. We separate the phonon field, which for a 
solid phases is a two component vector and a scalar for 
a smectic, into the high and small wavevector modes 



u(r) 



u 



; (r) 



u 



(4.1) 



and integrate perturbatively in Uk the high wavevector 
part u > (r) , with nonvanishing Fourier components inside 
a thin momentum shell 



ke~ l < |q| < A 



(4.2) 



We then rescale the lengths and long wavelength part of 
the fields with 



so as to restore the ultraviolet cutoff back to A = 2n/a. 
Because the pinning potential nonlinearity is a periodic 
function, it is convenient (but not necessary) to take the 
arbitrary field dimension to be 



= 0, 



(4.4) 



thereby preserving the period a — 2?t /A under the renor- 
malization group transformation. |44|] Under this trans- 
formation the resulting effective Hamiltonian, H = H e \ + 
i?K, can be restored into its original form with effective 
^-dependent elastic and Uk couplings. 

For the periodic pinning potential coupling Uk, we find 
in a standard way [m33] 



(4.5) 



U K {1) = U K e 2l ^ K »- 



where (w^)> is to be computed with the elastic Hamil- 
tonian appropriate to the phase being analyzed, keeping 
only modes within an infinitesimal momentum shell near 
the zone boundary A. Hence the nature of the pinning 
by the substrate potential and the transition tempera- 
ture obviously depend on the degree of the translational 
order in the system, i.e. whether the phase is a solid, or 
a smectic. 



1. Locked Floating Solid to Floating Solid Transition 



To determine the critical temperature for the LFS-to- 
FS transition, we compute the (u 2 } > a verag e using the 
anisotropic elastic Hamiltonian H e \, Eq.2.25, describing 



the 2d solid phase in the presence of a Id periodic po- 
j tential. Rewriting H e \ in terms of Fourier transformed 
phonon fields u(q), we find 



d 2 q 



(B xx q 2 x + B yx q 2 y ) \u x (q 



(K xy ql + K yy q 2 y ) |u y (q)| 2 + Sq x q y u x {q)u y {-q) 



(4.6) 



where, 



^XX 1 

/i-a + 7 , 

yy = ^yy ' 
K xy = // + a + 7 , 



B x 
K 



6 — fi + A 



x y 



■ 7 , 



(4.7a) 
(4.7b) 
(4.7c) 
(4.7d) 
(4.7e) 



which, after a simple Gaussian integration leads to 

k B T 



V> 



K, 



yy% 



K-xylx 



PqW y 
B^ql+Byxql 



where we have introduced the shorthand notation f = 

J y (2tt) 2 f° r ^ ne integral over the momentum shell. In 
the dilute limit and neglecting effects of the periodic po- 
tential on the elastic coefficients this reduces to 



«)> 



k B T 
2irJI 



(4.9) 



with ~p = 2[i(2[i + A)/(3/i + A). In order to compute 
( u y }> m general we use an elliptical (volume conserv- 
(4.8) ing) momentum shell defined by major and minor axes 

A X = kyfKyy/K X y SJld Ay = k^K X yjKyy. We find 
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k B Tci 



V> 



(4.10) 



27T^/ KyyK X y 

where we defined a dimensionless number c\ given by 



ci 



and 



2tt 



(W 



a x + {a y - a x ) 



o 2n a x + (a y - o z - a^y) sin 2 + a x?/ sin 9 ' 

(4.11) 



K. 

_ By 



By X 



yy 
5 2 



K K 

yy x y 



(4.12a) 
(4.12b) 

(4.12c) 
(4.12d) 



Upon combining Eq. 4.10| with Eq.4.5, we find the 
eigenvalue of the substrate potential to be 



k B Td 



o A, U _t Cl 

X P = 2-K 2 V 



(4.13) 



4TTy/KyyK X y 

which after setting 

A P (T p s) - (4.14) 
gives us the depinning transition temperature T p s 



8tt 



kBT pS = — yj KyyK X y ( — 



(4.15) 



which separates the LFS and FS phases. In the dilute 
limit the transition temperature reduces to 



(4.16) 



fcB^pS = 8tt/i — 

Z7T 



2. Locked Smectic to Floating Smectic Transition 



As discussed in Sec. Ill B, at low colloidal densities our 



system can exhibit LSm and FSm phases, and therefore 
undergo a phase transition between them in the roughen- 
ing universality class. Analogously to the LFS-FS tran- 
sition analyzed above, we can calculate the pinning tem- 
perature for the LSm-FSm transition by computing the 
(w 2 )> that goes into Eq j4.5| and finding the temperature 
at which this R G eig envalue vanishes. Using the Hamil- 
tonian Hs m , Eq. 3.35| , appropriate for the smectic phases 
and computing to zeroth order in the pinning potential 
Uk, we find 



V> 



k B T 



Byyl'y 



27T^/ ByyB X y 



(4.17) 
(4.18) 



where for convenience we again used an elliptical momen- 
tum shell with axes hyj B yy / B xy and K^J B xy /B yy . 

After combining this result with Eq.4.5 we find that 
translational pinning by the periodic potential is rele- 
vant in a floating smectic phase for T < T p g m , with T p g m 
given by 



k B T p sm — Bn^/ ByyB X y 



2n 



(4.19) 



As disc ussed in more detail in Sec.fv|, the clastic mod- 
uli in Eqs.l.lf; and 4. IE are functions of the strength of 



the pinning potential Uk, which in turn is proportional 
to the input laser intensity Ii n - Hence the resulting func- 
tions Tps(Iin) and T p s m (^in) in principle determine the 
LFS-FS and LSm-FSm phase boundaries displayed in 
Fig.|| for colloidal densities commensurate with the Id 
periodic potential. 



B. Dislocation Unbinding Transitions 



In the analysis of the preceding Sec.lVA, where we 
studied a thermal depinning transition within the solid 
phase, we implicitly assumed that the dislocations that 
distinguish the 2d solid and the smectic phases from 
the higher temperature disordered phases remain bound. 
Hen ce , these calculations for the pinning transition and 
Eqs. 4.15 , LIE remain valid only if they fall below the cor- 
responding dislocation unbinding melting transition tem- 
peratures, which we now compute. 



1. Locked Floating Solid to Locked Smectic Transition 

I t is e asy to see from the effective Hamiltonian -Hlfs, 
Eq.3.20, that the most striking consequence of the Id 
periodic potential is that it leads to the LFS phase, in 



which, the phonon degree of freedom, u y , correspond- 
ing to displacements transverse to the potential troughs 
acquires a "mass", Eq. B.1E, and as a consequence are ef- 
fectively suppressed. Therefore, this phonon mode can 
be safely integrated out, leaving an effective anisotropic 
2d XY Hamiltonian, with temperature and potential 
strength dependent effective elastic constants, 

H LFS = i J d 2 r[B yx {d v u x ) 2 + B xx {d x u x f) , (4.20) 

that describes a locked floating solid at sca les long er than 
the correlation length £ introduced in Sec. Ill AS . 

Melting of the LFS can be understood in terms of 
dislocation unbinding. However, in contrast to melting 
in the absence of an external (e.g., substrate or laser) 
potential, || here only the so-called type I dislocation 
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pairs (in the notation of R ef. pq) with Burgers vectors 
±t>i = ±bf{X (see Sec. |ll B| ) aligned parallel to the trough 
direction (which we continue to take along the x-axis) 
can be thermally unbound. In the presence of a periodic 
potential, oppositely charged dislocations, with Burgers 
vectors not satisfying the above condition (type II dislo- 
cations) are bound by a potential which grows linearly 
with the separation and therefore cannot thermally un- 
bind. This discussion is consistent with the mapping onto 
scalar Coulomb gas Hamiltonian, expected to describe 
logarithmically bound type I dislocations, embodied in 



the 2d anisotropic XY model Hamiltonian, Eq.4.20 



Away from the dislocation core, for a commensurate 
orientation defined by the shortest direct lattice vector 
pointing parallel to the troughs, R3 = n\ei + ri2e2, la- 
belled by dire ct lattice Miller indices n\ and n2 defined 
by Eq. 2.6-2.8, the displacement vector u for the active 



type I dislocation is given by 



x — tan 

2tt 



' D 1 / 2 ' 



xB 



1/2 



with 



bf L =| R s |= a J 7ii + n 2 + n i n 2 



(4.21) 



(4.22) 



Melting of the LFS via unbinding of these defects is iden- 
tical to the vortex unbinding transition of an anisotropic 
2d XY model. A standard calculation || leads to the 
prediction for the LFS melting temperature 



hi 



fcB?LFS-LSm = 77^" v B xx B yx 

on 



(4.23) 



and all other concomitant Kosterlitz-Thouless phe- 
nomenology This implies an exponential growths of the 
translational correlation length 0] 



6 



ae 



c/|T-T LFS 



1/2 



(4.24) 



with c a nonuniversal parameter and a universal ra- 
tio of the jump in the geometric mean of the shear 
and bulk moduli, Sj /I (T LFS _ LSm ) and B xx (^lfs— LSm) to 

TLFS-LSm- @ 

The resulting high temperature phase is the LSm jl6| , 
for low colloidal densities (i.e., high commensurability ra- 
tio p), and a modulated liquid for high densities (p < 1, 
see below), for which the smectic is indistinguishable 
from a liquid. Because of the unusually stro ng gr owth 
of the translational correlation length £ t , Eq. 4.24 , the 
phenomen ology o f the LSm-FSm transition that we stud- 
ied in Sec jlV A 2| will be modified for T -> T+ FS _ LSm by 
a long crossover from the crystal to smectic (or liquid) 
elasticity. 

It is important to note the distinction between this 
anisotropic 2d XY melting of a LFS into a LSm and 
an analogous type I melting mechanism of Ostlund 
and Halperin for melting of uniaxially anisotropic, but 



substrate-free 2d solids. j37j In the later case, thermal 
fluctuations destabilize the resulting 2d smectic by fur- 
ther unbinding type II dislocations, asymptotically con- 
verting it into a liquid. Here, because of the pinning 
potential, type II dislocations (e.g. ±b2,3 for p = (p, 0)) 
remain bound by a linear potential. The resulting LSm 
phase is therefore distinct from the (orientationally or- 
dered) modulated liquid (in which type II dislocations 
are also unbound), separated from it by a thermodynam- 
ically sharp phase transition. 



2. Floating Solid to Floating Smectic Transition 

A floating solid can melt continuously via unbinding 
of the type I dislocations. However, in contrast to the 
similar melting of a locked floating solid, here the dis- 
location unbinding in the displacement u x proceeds in 
the presence of another spectator massless phonon mode 



which is coupled to it. Consequently, as we will 
show below, this transition is a nontrivial extension of 
the Kosterlitz-Thouless theory and to our knowledge is, 
heretofore unexplored. Once these type I dislocations un- 
bind the most likely resulting phase is the floating smec- 
tic. @ 

The phenomenology of the FS-FSm melting transition 
can be most easily analyzed by the following steps. We (i) 
introduce dislocation degr ees of freedom into the elastic 



Hamiltonian H c \, Eq.2.25, (ii) perform a duality trans- 
formation to convert the resulting Coulomb gas Hamil- 
tonian into a modified Sine-Gordon model, and (iii) com- 
pute the dislocation unbinding temperature by analyzing 
the resulting dual model. 

To execute these standard steps, it is convenient to first 
perform the following rescalings of spatial coordinates: 



(4.25a) 
(4.25b) 



X - > X{B XX I By X ^j / 5 

y -> y{B yx /B xx )V 4 , 

which leads to the Hamiltonian 

Hfs = J d 2 r^K x (Vu x ) 2 + c x (d x ii y ) 2 + c y (dyUy) 2 

+ 2\ xy (d x u x )(d y Uy) + 2{n - "f){d x Uy){d y u x )} (4.26) 

where we dropped the prime on the rescaled coordinates 
and defined clastic constants 



-K X = \J B XX By X , 

c-x — -K x K X y / B xx , 

Cy — K X Kyyf By X . 



(4.27a) 
(4.27b) 
(4.27c) 



Because in the presence of dislocations the displacement 
field u x is a multivalued function, it is essential to dis- 
tinguish the last two terms in Eq, 4.26| . In contrast to 
conventional clastic theory, where dislocations are bound 
and u x is a well-defined function, here these terms can- 
not be transformed into each other by an integration 
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by parts. Keeping track of this distinction ensures the 
proper form for the elastic constants of the resulting 
smectic phase. 

In this new rescaled coordinate system, a type I dis- 
location located at the origin, with a Burgers vector 
b = 6,-jx, can be represented by a displacement field 



x — tan 
2tt 



(4.28) 



However, in contrast to the analysis of the melting of the 
LS phase above, in the presence of a finite d y u y deforma- 
tion, the form of type I dislocation given in Eq.4.28 does 
not correspond to a relaxed u x displacement which mini- 
mizes the energy. Consequently, we expect (see Eq. 1.34 ) 
a bilinear coupling between the dislocation density and 
the u y distortion. For a finite density of dislocations, we 
define a singular strain v s = VmJ due to a dislocation 
density b(r), with the standard relation 



V x v s 



E 



r), 

b H n rt S (2 \r- 



= e z b R n(r) 



(4.29) 
(4.30) 

(4.31) 



where the {n ri } are integer dislocation charges. A gen- 
eral solution to the above equation is given (in Fourier 
space) by 



v s (q) 



q 2 



6(q) + iqx(q) 



(4.32) 



where x(q) is an arbitrary, single- valued function, which 
for convenience and without loss of generality we can set 
to zero. After expressing the gradient of the total dis- 
placement field u t in terms of the dislocation part v s 
and a single valued phonon field u 



VUy , 



(4.33a) 
(4.33b) 



and inserting it into Hps , we obtain a Hamiltonian that 
includes both the elastic and dislocation degrees of free- 
dom 



H 



FSd 




X xyll + (7 - tilt 



H FS [u] . 



n(q)u y (-q) 
(4.34) 



After putting the system on the lattice, going to the 
grand canonical ensemble for dislocations, and adding 
the dislocation core energy E c to account for the energy 
coming from short length scales, (not included in above 
analysis), the total partition function is given by 



Z = [ [da] e~ BmA ~^ 

{n r } 



(4.35) 



In above, for convenience we chose to measure all the 
energies in units of k^T. 



To analyze the dislocation unbinding transition, it is 
convenient to perform a duality transformation on the 
above Hamiltonian i/FSd- |HJI] ^° ^° ^is we introduce 
an auxiliary Gaussian field (f> to decouple the Coulomb in- 
teraction between dislocations and use the Poisson sum- 
mation formula to perform the summation over the set 
of lattice integers {n r }, obtaining 



Z = J [d(j>][da]e- Hd 



(4.36) 



where, 
H A = 



VvMt + M)]} + H FS [u] . 

(4.37) 



To obtain H^, Eq. 4.37, above, we defined a field 0(r), 
whose Fourier transform is given by 



q- 



(A^ + (7-M)9x)%(q) . (4-38) 



and used Vv{4>) to denote the well-known 27r-periodic 
Villain potential defined by 



E 



n— — oo 



-E c n z +in4 



(4.39) 



At low fugacity (large core energy) , this potential reduces 
to a cosine function, leading to 

H d = [ d 2 r\ %H V</»| 2 - g cos [6 S (0 + %B)\ \ + H FS [a] , 



(4.40) 

with g = 2e- Ec . 

Now the dislocation unbinding transition in the origi- 
nal model of the floating solid is determined by the van- 
ishing of the RG eigenvalue of g(£) cosine nonlinearity in 
this dual model, defined by 

g{t) = ge^W , (4.41) 

where rj g is determined by 

Vg i = bl[{4> 2 ) > -{e 2 ) > ] , (4.42) 

with the right hand side easily computed fro m the 



quadratic part of the dual Hamiltonian H^, Eq. 4.37 
Specifically, 



= K r 



d 2 q 1 



2tt 



(4.43) 
(4.44) 



and 
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where, 



/ d 2 q {K v q; + (7 ~ M)gx) , ,,■_> 
(2tt) 2 g 4 



/ "TT^ , ".(q'l- (4.45a) 

(4.45b) 



d 2 q \ 2 xy q% + 2\ xy (~f - (i)qlql + (7 - ^) 2 g 4 



(2tt) 2 g 4 [ Cxg 2 + Cy(? 2 _ £2^2/(^2)] 



\\yCl + \ x y{j ~ fl)C 3 + (7 ~ M) 2 C4 . .„ . 

— - -^=^= I , (4.45c) 

111 KyyK X y 



f 2 * d9 a 2 sin 4 6> 

c 2 = / , (4.46a) 

Jo 2?r [a* + (a y - a K ) sin 2 0] [a^ + (a y — a x — a xy ) sin 2 + a xy sin 4 0] 

/" 27r d6 2a 2 cos 2 (9 sin 2 6 

c 3 = — , (4.46b) 

Jo 27r [a x + (ay - a x ) sin 2 6] [a x + (a y - a x - a xy ) sin 2 9 + a xy sin 4 6] 

f 2n de a 2 cos 4 

c 4 = / — T „ ^ —-2 ;— j — r- , (4.46c) 

J 27r \a x + (a y — a x ) sin 9\ \a x + [fl y — a x — a xy ) sin 9 + a xy sin 9\ 

Upon combining these results, we find that a floating solid melts into a floating smectic at 

8t V ^K yy K xy I 



wh ich r educes to the melting temperature Tlfs-lsih, 
Eq.4.23 of the LFS in the limit K xy , K yy — > 00, in which 
the spectator phonon u v mode is frozen out. Not sur- 
prisingly, we find that the extra u v fluctuations of the FS 
always suppress the melting temperature of the FS rela- 
tive to that of the LFS, i.e., for all range of parameters, 

^FS-FSm < ?LFS-LSm- 

We now demonstrate that once type I dislocations 
unbind, the resulting Hamiltonian is that of a floating 
sm ectic, described by the Hamiltonian Hpg mi given in 
Eq. 3.35. To see this return to the Hamiltonian -ffFScb 



Eq. |4.34| , and note that once dislocations unbind and 
therefore appear in large densities, the discrete disloca- 
tion field n r can, to a good approximation, be treated 
as a continuous density n(r). Within this Debye-Huckcl 
approximation the dislocation degrees of freedom can be 
easily integrated out of the partition function Eq.4.35 
by replacing the summation over n r in Z by an inte- 
gration. Simple Gaussian integrations over dislocation 
density n(r) and the single valued field u x then lead, in 
the long wavelength limit to an effective floating smectic 
Hamiltonian 



Hsm — 77 



d 2 , 



[n{d 2 x u y ) 2 + B xy {d x u y f + B yy (d y uy) 2 } , 



(4.48) 

where we have resto red t he original scaling of the spa- 
tial coordinates, Eqs. 4.25 and derived the effective elastic 
constants for the resulting FSm phase 



(M-7) 2 
(// + 7 - a) 2 



2E C 



(4.49a) 



_ 4/X7 - a 2 

13 X y — 

jtz + 7 - a 

\2 

R — X xy 

D yy *yy \ 

a x X 



(4.49b) 
(4.49c) 



We note that B xy vanishes as 7, a — > 0, as it must in 
this rotationally invariant limit, in which one must re- 
cover the rotationally invariant 2d liquid crystal smectic 
elasticity. j|| 

Another equivalent but considerably more straightfor- 
ward way to obtain the smectic Hamiltonian is to note 
that in t he pre sence of unbound type I dislocations the 
Vu* j Eq. 4.33a contains both the longitudinal and trans- 
verse components, and therefore, despite of its appear- 
ance is no longer a conservative vector constrained to be 
a gradient of a single-valued function. This observation 
allows us to incorporate unbou nd type I dislocations into 
the Hamiltonian H e \, Eq.2.25, by the replacement 



Vw^v, (4.50) 

with v an arbitrary 2d vector field. Under this substitu- 
tion H c \, Eq.2.25 transforms into 



tfFSd = j d 2 r{^d X Uy+Vy) 2 + ^V 2 x + ^f{d V Uy) 2 

+ KyV X dyU y + ^[{d X Uy) 2 ~ V^} + ~(d X <Uy ~ Vy)^ . 

(4.51) 

After performing a simple Gaussian integration over the 
two independent components of v, we immediately ob- 
tain a Hamiltonian for the floating smectic, which in the 
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long wavelength limit agrees in form and with the expres- 



sions for the elastic constants B T 



Eqs.4.48 and 4.49 



and B yy , obtained in 



V. SHAPE OF THE MELTING CURVE 
A. Strong pinning limit and reentrant melting 

One of the most interesting observations in the col- 
loidal experiments by Wei etal. ]2^ j, which in fact stim- 
ulated our interest in this problem, is the light-induced 
reentrant melting. As we shall explicitly demonstrate, 
this melting reentrance is a generic consequence of short- 
ranged screened colloidal interactions and thermal fluctu- 
ations, and hence should be prevalent in such 2d systems. 

To demonstrate the reentrance as a function of laser 
intensity, we study the shape of the melting curves for 
the LFS-ML, LFS-LSm and FS-FSm transitions, which 
we generally denote by T m (£/&-)■ The common fea- 
ture of these transitions is that they are all driven by 
the unbinding of type I dislocations, with T m (J7_R-) (see 
Eqs. 4.23 4.47 ) at least in part determined by the renor- 
malized values of the bulk modulus B xx for compression 
along the troughs and the corresponding shear modulus 
B yx . Our goal then is to detemine how these moduli 
depend on the potential amplitude Uk- 

We first note that these melting boundaries T m ([//f) 
are constrained by their limiting values 



T m (0) 



b* MM + A) 
4tt 2/i 
bl 



A 



T m (oo) = -VM2M + A) , 

S7T 



(5.1a) 
(5.1b) 



where T m (0) is the well-known result in the absence 
of an external potential pHq]. In the opposite li mit of 



infinite potential strength, T m (oo) is given by Eq. 4.23 
with B XX (U K 



oo) ps 2/i + A and B vx (Uk — » oo) 



Eq 



3.2C 



PS fj,. 
with 



These results follow from comparing Hlfs 
Ho, Eq. 2.4, after freezing out the u y degree of freedom 
(u y = 0) in Hq, as is appropriate in this Uk — > oo limit. 
Although in general there is no universal relation between 
T m (0) and T m (oo), in a dilute colloidal limit, relevant to 
the experiments of Wei etal. |B2|, t he two Lame coeffi- 
cients are equal, fi f» A, and Eqs.|5.1| reduce to 



2T(0) = M 



67T ' 



T^ 1 (oo) = y/3 fi 



1.3 7^(0) 



(5.2a) 

(5.2b) 
(5.2c) 



One might have thought that the melting temperature 
simply increases monotonically with Uk from T m (0) to 
T m (oo). However, as we will now show explicitly, the 
ity-mode thermal fluctuations, enhanced as the periodic 



potential is lowered from infinity, generically increase the 
melting temperature for na » 1. Consequently, the 
melting curve, T m ([/j<-), must have a maximum in this 
limit, implying reentrant melting for a band of tempera- 
tures as a function of the potential amplitude. 

The origin of the reentrance effect can be understood 
on a heursitic level as follows. Clearly, at small Uk, 
we expect that the increase in the strength of the pe- 
riodic potential suppresses thermal fluctuations in u y , 
thereby lowering the entropy of the liquid (or the smectic) 
state, and therefore making freezing into a lattice free- 
energetically less costly. This naturally leads to an in- 
crease of T^Uk) with Uk at low laser intensities. How- 
ever, for potential strengths Uk S> fcgT, this entropic 
contribution to the free energy becomes unimportant. In 
this large Uk limit, the behavior of T m (UK) is dominated 
by a different mechanism having to do with the reduction 
of the elastic constants with increasing Uk and decreas- 
ing temperature. To see this, note that the effective shear 



mod ulus B vx (Uk) which enters T m (UK) (see Eqs.4.23 
and 4.47 ), is determined by the screened Coulomb in- 
teraction, V(r) — Vq exp(— Kr)/r , between colloidal par- 
ticles in neighboring troughs. In order to find an effective 
shear modulus for the u^-modes, one needs to integrate 
out the massive modes corresponding to displacements 
perpendicular to the troughs of the laser potential. This 
will be the route taken further below. Heuristically, one 
should get roughly the same result by assuming that the 
dominant effect comes from the shear modulus B xy and 



simply averaging the potential over the massive u y 
grees of freedom, which yields 



de- 



B vx (Uk) 



-K|r n+ i- 



(5.3) 



where r„ and r n+ i are positions of nearest neighbor col- 
loidal particles belonging to the n-th and n + 1-st Bragg 
planes, running parallel to the laser potential troughs. 
This gives to lowest harmonic order in the fluctuations 



B vx (Uk) 



- na— k[u v (n+l) — u y (n)] \ 



e -rea e K (u y ) 

B yx (oc)e ck ° T / u « 



(5.4) 



with c a dimensionless number of order 1. Such a thermal 
enhancement of the effective shear modulus B vx (Uk), 
which decreases as thermal fluctuations in u y are sup- 
pressed by increasing Uk, is easy to understand: Even 
though, in the presence of u y fluctuations colloidal par- 
ticles in neighboring troughs spend as much time closer 
together as further apart, because of the concave form 
of the interaction potential the enhancement of the effec- 
tive shear modulus is larger from particles being closer 
together than the corresponding suppression when they 
are further apart. 

The above simple physical argument for reentrance is 
supported by detailed microscopic lattice calculations, in 
which we compute both the effective shear B vx (Uk) and 
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bulk B xx {Uk) moduli. To do this we start with a micro- 
scopic model with a screened repulsive Coulomb interac- 
tion V(r) — Vq exp(— nr)/r, where the screening length 
k _1 is typically much smaller than a and Vq depends 
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on the dielectric constant, k and the sphere radius p2[ . 
Upon integrating out the u y -modes using the screened 
Coulomb potential to leading order in k^T/UK the cal- 
culation in Appendix^ gives (for orientation p — (1, 0)), 



B yx (U K ) & B yx {oo) il + 
B XX (U K ) S3 B xx (oo) \ 1 + 
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where v — Vbe _Ka /fc^T ', B yx (oo) — ^vksTn 2 and B xx (oo) = 3^. Lowering the potential strength {Jr- always increases 
the shear modulus, whereas the behavior of the compressional modulus depends on the magnitude of v and na. When 
combined with E n. (j4~23| ), these expressions imply that the melting temperature increases with decreasing Ujc for 
na > 5.6 (in Ref. Ea/tow 10), 
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thus implying reentrant melting for a band of tempera- 
tures as a function of potential strength observed in ex- 
periments and illustrated in Figs .||4] and |[ Clearly given 
the dependence of the Tps-FSm on the elastic moduli, 
Eq. 4.47, we expect FS-FSm transition to display reen- 



trance, although quantitative predictions of the size of 
the reentrant band are much more difficult. 



In obtaining Eq. 5/7 we have clearly ignored addi- 
tional renormalization of the effective elastic constants 
by phonon nonlincarities and by bound dislocation pairs, 
which need to be taken into account for more precise esti- 
mate of the phase boundary. Based on general structure 
of Kosterlitz-Thouless-like RG flows, the latter renormal- 
izations generically reduce the elastic moduli and there- 
fore drive the melting temperature down. Since u y modes 
fluctuations and therefore the renormalizations that they 
induce are suppressed by the increasing periodic poten- 
tial, we expect that T m ([/ff) experiences larger reduction 
at small Uk than at large Uk ■ The known values for the 
the potential-free 2d melting and the 2d XY model down- 
ward renormalization constrain the extreme Uk = and 
Uk —* oo ends of the melting curve. Furthermore, since 
thermal downward renormalization of elastic constants is 
obviously enhanced with increasing temperature, we ex- 
pect the suppression of the melting temperature due to 
these effects to be most pronounced near the maximum in 
T-ox{Uk)- Clearly, such a [/^--dependent downward renor- 
malization of the elastic constants will generically tend 
to reduce the range of temperatures over which there is 
laser-induced reentrant melting. However, these effects 
are small 47 and we therefore expect reentrant melting 
to persist even in their presence. 



B. Weak pinning: universal shape of the melting 
curve at small potential strength 

In addition to a maximum displayed by the melt- 
ing curve as a function of laser intensity, we also find 
that the shape of the melting temperature is universal 
in the limit of a vanishing periodic potential strength 
Uk- This can be seen most easily from the RG scal- 
ing theory applied to the potential- free critical point. 
More specifically, consider the behavior of the transla- 
tional correlation length £(t, Uk) above the melting tran- 
sition as a function of Uk and the reduced temperature 
t = [T — T m (U K = 0)}/T m (U K = 0). The power of 
the renormalization group transformation is that it al- 
lows us to relate a difficult calculation very close to the 
transition, where fluctuations are large and perturbation 
theory is divergent, to a calculation outside of the critical 
region, where perturbation theory is convergent. Apply- 
ing this idea to the computation of £(i, Uk) we find 



t{t,U K ) ^b^{t{b*),U K b x , K ) , 
= e c / t77 e(l,C/Ke cAA " /tT7 ) 



(5.8a) 
(5.8b) 



where we have chosen the RG rescaling parameter 6* such 
that the rescaled reduced temperature i (6*), given by the 
RG flow equations of Halperin and Nelson || is of order 
unity 



t(M = i 



(5.9) 



Aa' = 2 — ?7a/2 is the renormalization group eigenvalue of 
the Id periodic potential Uk At the primary potential- 
free fixed point with Uk — 0, we recover the well- 
known |j exponential growth of the correlation length 
£(t, Uk = 0) with the exponent V given by 



0.36963 



(5.10) 



where an overbar denotes critical exponents at this fixed 
point. 
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The primary critical behavior is unstable for arbitrarily 
small Uk- Hence sufficiently close to the melting temper- 
ature T m (0), the periodic potential always becomes im- 
portant. This is the case even for the melting of the FS, 
where it leads to a marginal crossover from a fixed line 
of isotropic rotationally invariant elastic ity t o the fixed 
line characterizing elasticity given by Eq.2.25, where the 
rotational symmetry is explicitly broken by Uk- There, 
despite the fact that the periodic potential is irrelevant 
for the translational order parameter, it is always impor- 
tant for the orientational degrees of freedom, since (see 
Sec. IIC) it explicitly breaks orientational symmetry. 

In locked phases, it is clear from Eq. 5.8b, that for a 
given small Uk, the effects of this weak periodic poten- 
tial will be felt at T m (U K ) > T m (0), such that t he Uk - 
dependent argument on the right hand side of Eq.5.8t is 
large, i.e. grows beyond order T m , 



Uk ~ kBT m e 



(5.11) 



This then predicts a universal cusp for the melting curve 
T m (UK) in the limit Uk — ► in any phase in which 
Xk > 0, i.e., the periodic potential is relevant and the 
phase is locked and T m (UK ) given by 

T m (U K ) ~ T m (0)[l + XK[Hk B T m /U K )]- 1/17 ] , (5.12) 

as depicted in Figs.0,0 and 0. For floating phases, such 
as the FS and FSm, where the periodic potential is irrel- 
evant (in the RG sense), we expect the convergent per- 
turbation theory in Uk to lead to a melting temperature 
(£/#■) that instead grows linearly with Uk- 



VI. RESPONSE OF THE TRANSLATION AND 
HEXATIC ORDER PARAMETER TO AN 
EXTERNAL POTENTIAL 

In this section we use a renormalization group scaling 
analysis to determine the response of the translational 
order parameter Mk = (pa) an d the bond orientational 



order parameter = (e 



6i0(r)\ 



to the amplitude Uk of 



the external laser potential. In the absence of an external 
potential, Uk = 0, there are only algebraic peaks in the 
static structure function of the crystalline phase and the 
translational order parameter Mk = (pk) vanishes like 
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as the system size L — > oo, where 
kuT 3/i + A 



4tt /z(2/z + A) 
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(6.1) 



(6.2) 



is the critical exponent of the potential-free case §. For 
small values of the external potential Uk we can use 
standard crossover scaling analysis to determine how the 
translational order parameter depends on the amplitude 



of the laser potential. We start from the scaling behavior 
of the free energy density under a renormalization group 
transformation 



f(U K ,T) 
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f(e x « e U K ,T(l)) , 



(6.3) 



where Ak is the renormalization group eigenvalue for the 
periodic potential, and T{tj is the renormalized temper- 
ature which characterizes the crystalline phase. Since in 
the free energy density the laser potential Uk couples 
linearly to pc we have 
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K 



dU 
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■f(U K ,T) 



(6.4) 



and dimensional analysis tells us that the exponent of 
the correlation function (pkCOpkW) ~ r~ r ' K is related 
to A K by 



A 



K 



(6.5) 



a result consistent with standard perturbative calcula- 
tion of Ak- Hence we get the following scaling relation 
for the translational order parameter 

M K (U K ,T) = e-i^ e M K (e^-^ e U K ,T(e)) , (6.6) 

where we expect T(i) to approach a finite value as 
t — ► oo. Upon choosing I = £„ such that e( 2_77K ' 2 "* Uk — 
/ia 2 , i.e., is comparable to the elastic energy for deforma- 
tion at the lattice cutoff a, we obtain 



Mk(U k ,T) H U k 



(6.7) 



For ?7 K > 2, Mk vanishes linearly with Uk, with a sin- 
gular correction. In contrast, Mk should always vanish 
linearly with Uk in the liquid and hexatic |li"| , |l2| ] phases 
of the unperturbed colloid. 

The laser potential will also induce long-range bond 
orientational order in tpe — (e 610 ^) fig . Along similar 
lines as above, one can show that the bond order param- 
eter ipe vanishes linearly with Uk in the liquid, vanishes 
like a power of Uk in the hexatic phase 



A H U K I 



6v 6 /(4-V 6 ) 



(6.8) 



where rj e is the exponent describing the algebraic de- 
cay of bond order, and approaches a nonzero constant 
as Uk — > in the solid phase, pq] 



VII. DISCUSSION AND EXPERIMENTAL 
IMPLICATIONS 

A. Melting temperatures and critical 
commensurability ratios in the dilute limit 

One of the interesting predictions of our work is that 
the LFS-ML, LFS-LSm and FS-FSm transition are all 
mediated by the unbinding of type-I dislocations with 
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Burgers vectors parallel to the troughs of the external 
potential, b = 6^x. Consequently, depending on the 
choice of relative orientation, the periodic potential can 
be used to supress the unbinding of a set of dislocations 
that would otherwise unbind in a "substrate" -free ex- 
periments. For example, in the dual-primary orientation 
shown in Fig.|l^, all six fundamental Burgers vectors are 
confined by a linear potential and therefore cannot un- 
bind entropically. It is therefore the unbinding of non- 
fundamental dislocations with Burgers vector of charge 
V3a, illustrated in Fig.|l7] that will control the melting 
transition. 




FIG. 17. Triangular lattice with lattice constant a subject 
to a periodic potential (maxima indicated by dashed lines) 
for a dual-primary orientation with pd = a', where a' = a/2 
is the Bragg plane spacing and the commensurability ratio is 
p — 1. Also shown is the low energy dislocation with Burgers 
vector b parallel to the corrugation of the potential. 



In general, the magnitude 6^ = a 2 (n\ 



nin 2 ) 



of the lowest energy Burgers vector and hence the melt- 
ing temperature T m oc 6| depends strongly on the rela- 
tive orientation between the 2d solid and the laser po- 
tential, e.g., for n% = 1 and n\ = 0, 1, 2, 3 one finds 
b\ja 2 = 1, 3, 7, 13. In particular, if one keeps the mean 
particle spacing a (i.e., the density) and the potential 
strength fixed and reduces the spacing d between the 
laser troughs (by e.g., varying the angle between the 
two interfering laser beams) d — > d/V3 such that one 
goes from preferred lattice orientation A (ua = (1,0)) to 
preferred lattice orientation B (ris = (1,1)) the melting 
temperature should increase by a factor of 3 (see vertical 
arrow in Fig. |l8|a). This appears to be consistent with 
preliminary data of Bechinger etal. E^j. They find that 
for exactly such a change in trough spacing the onset of 
light induced freezing at fixed temperature is shifted to 
smaller laser intensities also by roughly a factor of 3 (see 
horizontal arrow in Fig. |l8|a). 

More detailed experimental studies of T m {UK) for 
various commensurate orientations and trough spacings 
would clearly be desirable in order to systemmatically 
test our predictions for the orientation dependence of 
the melting transition temperature for the LIF. In per- 
forming such studies one must keep in mind considerable 



irreversibility effects that are expected to plague "zero- 
laser-field" cooled experiments. In order to avoid dealing 
with long equilibration times, one would need to warm 
up into the liquid state, change the laser potential period 
d and only then "field-cool" back into the solid. 
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FIG. 18. a) Schematic p — 1 phase diagram as a function 
of potential strength Uk and relative orientation between the 
laser potential and the 2d solid; a change in orientation from 
(1,0) to (1,1) is generated by keeping the particle density 
(and hence the mean particle spacing a) fixed and varying 
the distance d between the minima of the external potential, 
b) Schematic p — 1 phase diagram for lattice spacing d ver- 
sus potential strength Uk at fixed temperature T and fixed 
colloidal density; incommensurability effects are disregarded. 

Since trough spacing d (controlled by the angle be- 
tween the interfering laser beams) and laser intensity ap- 
pear to be convenient experimentally tunable parame- 
ters, it is valuable to derive the shape of the melting 
curve in the d-Ux plane (for a given temperature and a 
fixed density of colloidal particles). However, since an 
arbitrary value of d will in general not be commensurate 
with the spacing between a particular fixed set of Bragg 
planes, a detailed study of incommensurate potentials 
would need to be done in order to fully understand the 
behavior as a function of trough spacing d. We hope 
to discuss some of the ensuing physics in a forthcoming 
publication p7[ |. But, for the following we would like 
to restrict ourselves to values of d which are commensu- 
rate. Hence, strictly speaking, our results will be not be 
valid for a continuous set of layer spacings but only for a 
discrete commensurate subset of values. With this pre- 
caution in mind we expect the melting curve (for a given 
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temperature and particle density) in the d-Ux plane to 
have the shape illustrated in Fig. 18 b. We note that in 



the LIF regime the critical potential strength for melt- 
ing decreases with decreasing distance between the laser 
fringes, whereas in the LIM regime the critical poten- 
tial strength increases as the interference fringes become 
narrower. 

Let us now specialize to the dilute limit, ita> 1 , rel- 
evant to the experiments of Wei etal. p^ |. Then the two 
Lame coefficients (characterizing the continuum elastic 
theory of the hexagonal crystal in the absence of a laser 
potential) become equal, p ~ A, and the melting temper- 
ature for the LFS reduces to 
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(7.1) 
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in the limit of zero and infinite potential, respectively. 
For small values of the commensurability ratio, p < p c , 
the LFS melts into a modulated liquid or a locked smec- 
tic. If p > p c , a floating solid with two soft phonon modes 
can intervene between the LFS and a modula ted li quid 
or floating smectic phase. As discussed in Sec. IV A , the 
transition from the LFS into the intermediate FS phase is 
in a roughening universality class where the laser poten- 
tial becomes irrelevant. In the dilute limit (and neglect- 
ing effects of the periodic potential on the elastic coeffi- 
cients) the corresponding critical temperature is aproxi- 
mately 



(7.3) 



where d = a! jp with a! ~ \[Zaj (2y/nf + n| + n\n?) the 
distance between the Bragg planes parallel to th e tro ughs 
of the laser potential. Upon combining Eq. 7.3 with 
Eq. 7.1, the critical commensurability ratio reads 
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Note that only for the primary (p^ 11 = 3-\/3/2 w 3.7) and 
for the dual-primary orientation (p^ 1 = \/3/2 w 1.2) is 
this critical value larger than 1. For any other orienta- 
tion p c is less than 1 and hence we expect that there will 
always be an intervening floating solid pha se. A con- 
figuration with ft — (2,1) (see also Fig. [l(j) and hence 
p c = (3/7)-\/3/2 « 0.5 is likely to be within the range of 
parameters accessible to experiments with colloidal par- 
ticles. 

For p > p c there is a roughening transition from a 
locked floating solid into a uniaxially anisotropic float- 
ing solid described by flei, Eq. 2.25 , which subsequently 
melts (by unbinding of type I dislocations) into either a 
modulated liquid or a floating smectic. Since the melt- 
ing and the roughening transition for a locked smec- 
tic phase are respectively given by T m s m = i^Ba 2 and 



TpSm = -^Bd 2 , where B = B xy B yy and a is the smectic 
layer spacing, there exists a universal commensurability 
ratio p' c — 4 |4|] above which a floating smectic phase 
intervenes between a locked smectic or a floating solid 
and the modulated liquid. This universal value p' c = 4 
should be contrasted with the nonuniversal critical com- 
mensurability ratio p c for the existence of the floating 
solid phase, which depends on the relative magnitude of 
the elastic constants and strongly on the relative orien- 
tation between the colloidal lattice and the Id periodic 
potential. Current experiments find it difficult to access 
large commensurability ratios p. We hope that our the- 
oretical results will inspire experimentalists to overcome 
present obstacles and map out the rich phase diagram 
shown in Fig.||. 



B. Phase diagrams as a function of the Debye 
screening length 

Recent Monte-Carlo simulation studies of melting in 
the presence of a Id periodic external potential have 
explored the phase diagram in the parameter space of 
Ux/ksT and na with particle density and temperature 
fixed. 1 19[ 2(| Although one might question whether such 
simulations are in equilibrium with respect to disloca- 
tion climb (or even glide), it is important to tabulate the 
predictions of our defect-mediated melting theory in this 
parameter space in order to be able to compare with the 
results of these simulations. In addition, it also seems to 
be more feasable experimentally to map out the phase 
diagram as a function of potential strength and particle 
density. 

Adapting our results from Sec.^ we find the following 
behavior. Since the melting temperature is proportional 
to the elastic moduli, which in turn are proportional to 
the potential strength, for /ta > 1 we expect T m to dis- 
play the following dependence on the screening length 



T m oc (na) e Ka . As an immediate consequence one gets 
(in the dilute li mit) the following implicit equation (see 
also Eq. O and O) 



(«m - K m)« « 2 In (1.3 
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In particular this implies that the difference in the critical 
values of the inverse screening length at infinite and zero 
potential stren gth, and k^, is positive. In the limit 

<)o w 21nl.3 « 



.0 



K^a ^> 1, Eq. 7.5 reduces to (k° 

0.52. The full solution of Eq. 7J3| together with the 



asymptotic result is shown in Figjl9|. We find Eq. 7.5 
to be consistent with experimental results p9[ . It would 
be interesting to test experimentally the functional de- 
pendence of on K^, Fig.[l9|, predicted here. 

The results of Monte-Carlo simulations appear to dis- 
agree with experiments and with our predictions from 
the dislocation-mediated melting theory when compared 
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for large values of the potential strength. Whereas we 
find > k^, the simulations reported in Ref. Jll| show 
quite the opposite. More recent simulations from the 
same group [£(J seem to refute these earlier results and 
find in agreement with our theory — > 0. Their 
numerical value for (n^ — K^)a ~ 1.32 is, however, more 
than two times larger than our asymptotic prediction of 
0.52. However, because Eq. 7.5 neglects finite renormal- 
ization of elastic constants by dislocation dipoles and 
nonlinear elastic effects, our prediction is an estimate, 
only accurate upto unknown factors of order 1. 
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FIG. 19. Difference between K m a at infinite and zero po- 
tential strength as a function of K^a. The vertical dashed line 
gives the asymptotic value 2 In 1.3 ~ 0.52 for very large K^a. 

Next we discuss reentrance in t he Uk /k^,T-K, m a phase 
diagram. Upon rewriting Eq. 5.7 we find 



a(K m a) 



Uk = 
k B T T/T™{ Km a) - 1 



with 
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Hence, if K^a and n^a are both smaller than the crit- 
ical value 5.6 for the existence of reentrance, we expect 
(«; m a) _1 to be a monotonically decreasing function of the 
potential strength, as shown by the dashed line in Fig.|2C|. 
If K^a and n^a are both larger than the critical value 
5.6, we expect reentrant behavior such that with increas- 
ing potential strength (K m a) _1 first decreases and reaches 
a minimum (K™ m a) _1 < (k^ci) -1 before it approac hes 
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(k^ci) -1 as an inverse power of Uk according to Eq 
(see Fig.^o|). This reentrant behavior is consistent with 
results from experiments of the Konstanz group [ p2|j49| 
(see the dashed arrow in Figj2^, which describes a typical 
experimental path) . It is also similar to what one finds in 
simulations |l9[ at small values of the potential strength. 
However, there are significant differences. First of all, the 
type of transition is very different. Whereas we discuss a 
continuous dislocation mediated melting transition, sim- 
ulations appear to find a first-order transition. Second, 



as discussed above, the simulations show < which 
is opposite to what our theory predicts. In more recent 
simulations pp| K m a is found to increase monotonically 
with potential strength with no sign for reentrance. This 
is opposite to what was found in the earlier simulations 
by the same group lll| . 
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FIG. 20. Schematic p = 1 phase diagram as a function of 
potential strength Uk and inverse Debye screening length k. 
Solid and dashed curves represent the melting curves for val- 
ues of k being larger or smaller respectively than the critical 
value of Kcritffl ~ 5.6. 

In summary, we find that our theoretical results are 
consistent with recent experiments and raise strong 
doubts on the validity of the Monte-Carlo results to date 
on melting in a Id periodic potential. This latter failure 
of simulations is not completely surprising given difficul- 
ties of numerical methods on even larger systems to re- 
solved the nature of 2d melting even without an external 
potential fl2lf. 



C. Static structure factor and pair correlation 
function 



The quantity that is most directly observed in many 
experiments on colloidal systems and related simulations 
is the pair correlation function, defined by 
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where the double sum is over N particles but excludes 
the diagonal terms where i = j. It is related to the static 
structure factor by 

9{r) =^J S giq r ^ (q) ■ (7 - 9) 

Neglecting the smooth part of the structure factor and 
taking into account only the center column of Bragg 
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peaks and the two neighboring columns of quasi-Bragg 
peaks with r] GAl: (T m ) = 1/4 and r)G B ,A T m) = 1 > respec- 
tively, one finds for the pair correlation function, 

<?(r)-l= 5^C Gi cos(Git) 

; C Ga 008(0,4 • r) 

Ga 

-"^^C Gb cos(G b t), (7.10) 
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where nc a x are the exponents characte ristic for the LFS 
phase. Note also that according to Eq. 1.13 these expo- 
nents do only depend on the x— component of the recip- 
rocal lattice vector G a . The amplitudes Cg q are pro- 
portional to the amplitudes of the corresponding Bragg 
peaks, Gi = (2Tr/a')e y , and quasi-Bragg peaks, Ga and 
G B with G^x = G° x = 2tt/cl and G B . X = 2G° X (see 
Eqs. |3.32a| - ^33b| ). 

For r || x, i.e., looking parallel to the minima of the 
troughs the sum over the Bragg peaks yields a constant. 
This simply reflects the effect of the laser potential to in- 
duce a periodic modulation of the colloidal particle den- 
sity with a higher density in the minima of the troughs. 
Note that this trivially implies that the pair correlation 
function does not approach unity as x — > oo if g{x) is 
normalized with respect to the mean density. Since the 
amplitudes for the quasi-Bragg peaks decay as a power- 
law in the strength of the laser potential with an expo- 
nent proportional to rjQ y a reasonable approximation for 
the pair correlation function reads 



g(x) — 1 = const. + gACOs(G x x)x vg a, i 
+g B cos{2G x x)x~ riG z 



(7.11) 



The relative magnitude of the amplitudes gA and g B de- 
pends on the strength of the laser potential. Whereas 
g B is independent of Uk (note that the leading quasi 
Bragg-peak contributing to g B has G BtV = 0), gA van- 
ishes as a nontrivial T-dep endent power law in Uk for 
Uk /^o? -C 1 (see Eq. 3.32t| ) increasing the weight of the 



x ,g a,x -term with increasing potential strength. This 
prediction should be accessible to experimental verifica- 
tion. Note, that the dependence of the amplitude gA on 
the potential strength may lead to [/^-dependent effec- 
tive exponents when one tries to incorrectly fit the ex- 
perimental data by a single power law. For illustration 
Fig.|l| shows g(x) — 1 for a special case, where const. = 0, 
.9.4 = 9b = 1, r/G A ,* = ji VGb,* = 1 and au length are 
measured in units of a. Due to the superposition of the 
two harmonics with different power law amplitudes the 
minima are much broader than the maxima of the struc- 
ture factor, a feature which appears to be present in the 
data of Ref. & 




x 



FIG. 



21. Algebraic part of the static structure factor for 
r = x, i.e. looking parallel to the troughs of the laser poten- 
tial. There are two oscillating contributions, both of which 
decay algebraically to zero. 




y 

FIG. 22. Algebraic part of the static structure factor for 
r = y, i.e. looking perpendicular to the troughs of the laser 
potential. Unlike Fig.Ell algebraically decaying oscillations are 
superimposed on a periodic contribution (not shown) which 
does not decay. 

For r || y, i.e., looking perpendicular to the minima of 
the troughs we get, 

g(y) - 1 = const.' cos(2G° y y) + ^cos^y)*/" 7 ^.* 

+9 B y~ naB '° (7-12) 

with Gy = 2n/y3a. Hence on top of the periodic den- 
sity modulation due to the laser potential we have again 
the algebraic decay from the closest Bragg peaks. For 
illustration Figj2^ shows the algebraic part of the static 
structure factor f(x) = y^ 1 ^ 4 cos(2ny/^/3) + y^ 1 where 
we have again chosen the amplitudes to be equal and the 
n exponents equal to their values at the melting temper- 
ature, no A „ = 2, lG Sj = 1. If one would try to fit the 



31 



envelope of this function in the regime shown in the graph 
using a single power law one would find an exponent of 
1/2. Hence caution must be exercised in the analysis 
of the experimental data, and it is essential to take into 
account both leading and subleading quasi Bragg peaks. 
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APPENDIX A: VARIATIONAL THEORY OF THE 
2D MELTING TRANSITION IN THE PRESENCE 
OF A ID PERIODIC POTENTIAL 

In this appendix we study the freezing transition of 
the modulated liquid in the limit of a strong periodic po- 
tential. In such limit the colloidal particles are tightly 
confined to the troughs of the Id periodic potential and 
our system reduces to a weakly coupled array of Id col- 
loidal liquids. The low energy degrees of freedom of the 
resulting system are then well characterized by a scalar 
field u n {x) describing particle displacements along the 
ri-th trough and an effective Hamiltonian 



gcos [(f> n+ i(x) - (f>„(x))] 



(Al) 



where for simplicity of notation we have defined rescaled 
phonon field <p n (x) and elastic couplings B and g related 
to those defined in the Introduction through 



(j>n{x) 



2tt 



u n (x) , 



B = Kd 



g = fid 



— V 

2nd J 



(A2a) 
(A2b) 
(A2c) 



In the Introduction we have used simple qualitative ar- 
guments to estimate the colloidal freezing transition tem- 
perature. Here we would like to treat this model quan- 
titatively and in more detail. Unfortunately, however, as 



can be seen from a standard renormalization group anal- 
ysis weak coupling g is always irrelevant at long scales, 
with the effective coupling g(£ x ) vanishing at length scale 
S,x as 



(A3) 



Thermal fluctuations, which are especially strong in 
Id are responsible for this effective decoupling of 
the colloidal system into effectively independent one- 
dimensional liquids. This precludes a description of the 
freezing transition in weak (g) coupling starting from this 
model. There are two alternatives: One is to study of 
the melting transition from a complementary strong cou- 
pling, fully elastic model with topological defects (dis- 
locations), an approach which lends itself to a rigorous 
treatment that we undertake in the main part of the 
paper. Alternatively, an approximate, variational treat- 
ment of the model, Eq. |Al| is possible and will be pre- 
sented in this appendix. 

The idea behind a variational approach of a problem 
is that an approximate free energy 



(H - H„ 



F, 



(A4) 



is an upper-bound for the exact free energy F correspond- 
ing to the Hamiltonian H of interest and where H v any 
other (the so called variational) Hamiltonian, F v is the 
corresponding free energy and subscript v on the thermal 
average indicates that Boltzmann weight with Hamilto- 
nian H v is used. The advantage of the variational prin- 
ciple can be taken if the arbitrary variational Hamilto- 
nian H v is judiciously chosen to be simple enough, so 
that thermal averages can be calculated, but at the same 
time general enough so as to be able approximately cap- 
ture the physics of the full Hamiltonian H . 

Since, unfortunately, our abilities to compute func- 
tional integrals do not extend beyond Gaussians, we 
choose a quadratic form for H v 



H - = J2j dx 



F>x 
2 



dx 



B„ 



4>n) 



(A5) 



with B x and B y as the effective variational parameters, 
respectively related to the effective long wavelength bulk 
and shear moduli, latter given by 



/i = B y d 



2tt 



(A6) 



Simple Gaussian averages then lead to the variational 
free energy density f(B X) B y ) = F(B X , B y )/L x N y 
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/ = 



-(B - B x )kl - By{l - COS kyd) 



~k B TlogG v (k) 



G„(k) 



-ffexp 



— / (1 — cos kyd) G v (k.) 



(A7) 



where L x and N y (= L y /d) are respectively the length 
and the number of laser potential troughs (i.e., the 2d 
dimensions of our colloidal system) and G v (k) is the 
Fourier transform of the intra-trough displacement cor- 
relation function given by 



G„(k) = k B T [B x k 2 x + 2B y {\ - coskyd)}' 



(A8) 



To find the upper-bound of the free energy density /, 
we now minimize f(B x , B y ) over the variational parame- 
ters B x and B y . Conceptually simple but tedious calcu- 
lation gives 



B x — B 
B y( B -,9) = 9 e 



-k B T/[Kd(B y B) 1/2 } 



(A9a) 
(A9b) 



Equation A9b, which determines the behavior of B y 
and therefore the effective shear modulus \x as a function 
of temperature and inter-trough coupling g, illustrated 
in Fig.E3 is the main result of the variational calculation. 



Simple graphical analysis of Eq. A9b predicts 

B y {g) = 0, forg < g c , (AlOa) 
B v (g) w ge-^K^ 1 ' 2 , for g » g c , (AlOb) 



where, the critical value of the coupling g which sepa- 
rates the two solutions for B y is given by 



9c 



k^Te 
2tt 



1 

B 



(All) 








liquid g c crystal g 



FIG. 23. Shear modulus /i as a function of the inter-trough 
coupling g (at fixed temperature), showing a freezing transi- 
tion between /i = 2d liquid and a /i > 2d crystal, and a 
jump discontinuity at g c in the shear modulus. 



Combining this with Eq. AE , wc conclude that the tran 



sition between the two solutions in Eq. A10 represents the 
freezing of a zero shear modulus (/i = 0) 2d liquid into 
a finite shear modulus (// > 0) 2d solid. In terms of the 
she ar mo dulus /i and the bulk modulus K, defined by 
Eqs.A2c the corresponding melting transition tempera- 



ture is given by 



kBT m = — 1 

Z7T 



(A12) 



J 



a value that, upto factors of order 1, is consistent with 
the asymptotically exact prediction of our strong cou- 
pling (elastic model) analysis given in the main text. 



APPENDIX B: EFFECTIVE ELASTIC 
CONSTANTS FOR SCREENED REPULSIVE 
COULOMB POTENTIAL 



To calculate the effective clastic constants in the limit 
of large through potential, we start from a model with 
a pair potential given by a screened repulsive Coulomb 
potential V(r) — Voaexp(— nr)/r, where the screening 
length k -1 is typically much shorter that the mean par- 
ticle spacing a. The total potential energy is then given 
by 



i 



Rw 



(Bl) 



where due to the short range of the potential we can 
safely restrict summation to nearest neighbors, (1,1'). 
The distance between the colloidal particles numbered 
I and V can (for a perfect lattice) be decomposed into 
a distance between the equilibrium positions r; and the 
displacement vectors u/: 



Hw = ri- t v + U; - Up 
= r U i + u w 



(B2) 



In the following we restrict ourselves to the primary con- 
figurations and write the potential energy as sums over 
Bragg "planes" (i.e., rows of particles in d = 2) indexed 
by an integer r and particles within these rows indexed 
by^, 



[(a + Sui) 2 + Sh 



2lV2 



exp 



■Ac [(a + 8uiY 



Shf 



1/2' 



33 



1 



[(a/2 + Aw,) 2 + (d + A/i ; ) 2 ] 1/2 
1 



exp 



exp 



- K [(a/2 + A U/ ) 2 + (d + A/i/) 2 ] 1/2 
-K [(a/2 + Aw,) 2 + {d + Aht) 2 ] 1/2 ] } 



[(a/2 + A?20 2 + (rf + A/i i ) 2 ] 1/2 
where the relative intra- valley and inter- valley displacement fields are defined as follows (see Fig.p4|): 

Sui = u x (xi + a,yi) - u x {xi,yi) , 
Shi = u y (xi +a,yi) - u y {xi,yi) , 

and 

Aui = u x {xi + a/2, yi +d)- u x {x u yi) , 
Aui = -u x (x t - a/2,y t ) + u x (x h yi) , 
Ahi = u y (xi + a/2,yi + d) - u y (xi,yi) . 



(B3) 



(B4) 
(B5) 

(B6) 
(B7) 
(B8) 




1+1 



8uj 

FIG. 24. Sketch of two rows of a triangular lattice of colloidal perticles in a through potential illustrating two contributions 
to the effective potential energy. The sum over the lattice sites is done by summing along the valleys; there is one intra-valley 
nearest neighbor (wiggly line) and two inter-valley nearest neighbors one in the forward direction (Aw) and one in the backward 
direction (Aw)) (solid lines). 



In the strong pinning limit the laser potential can 
be expanded in powers of the phonon fields in the y- 
direction, 



where we have used pd = ^fia/2. In the following we 
shall (in order to simplify notation) measure all lengths 
in units of the mean lattice spacing. 



Uk 
k B T 



E 



277 u 
COS | — -til 



2 87T^ U K x - / U4 

P 3 k B T V ' " 



k B Tw V* ( — 



(B9) 



We proceed as follows: (i) first we expand all terms 
in the total potential energy <E> to quadratic order in the 
out-of valley displacement fields, (ii) integrate out the 
massive out-of-valley modes, and (iii) take the contin- 
uum limit. Note, that it is only step (i) which explicitly 
depends on the particular form of the pair potential. For 
simplicity, we will limit our derivation to the leading or- 
der in Vq/Uk and k B T/U K . 



Step (i) gives 



with 



$i[u,h] 
<f> 2 [u,h] 



f3<S>[u,h] =/3(#i +$2 + $s) 



v^2[-^(k + l)Shf + i( K 2 + 2k + 2)Suf + (3i(K)Sh 2 Sui + S^ShfSuf^ 



5k + 5)Ahf + -(k 2 - k - l)Auf + a 2 {n) ShiAui 



+f3 2 (n)AhfAu l + 7 2 (k)A/i ; Am 2 + <5 2 (k)A/i 2 Am 2 j 



(BIO) 



(Bll) 



(B12) 
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and $3[m, h] obtained from $2[«, h) by the replacement 
Aui — > Aw; . Here we have also introduced 



and 



ft («) = 2( k2 + 3k + 3 ) . 

<5i(k) = -|(« 3 + 5k + 12k + 12) 



\/3, 2 



Qf2 («) = -^-(k 2 + 3k + 3) 
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(3k 3 + 14k 2 + 33k + 33) 



72(«) = -^g( K3 + 2k2 + 3k + 3) 



(B13) 
(B14) 

(B15) 
(B16) 

(B17) 



S^k) = g|(3K 4 + 14k 3 + 55k 2 + 123k + 123) . (B18) 



The dimensionless ratio 



v = e 



(B19) 



measures the strength of the pair potential relative to a 
typical thermal energy. Next we integrate out the mas- 
sive phonon fields hi with a Boltzmann weight given by 
the external potential -ffic, 



cxp [-/3i/eff] = J[dh] exp 



, (B20) 



where J[dh] denotes an integration over the {hi}. We 
find 



i 



J(k 2 + 2k + 2) + -Ji(/c) - — o|(k) 
2 



+ (A« 2 + Au 2 ) [^(k 2 - k - 1) + ^ 2 (k)] } (B21) 

In the continuum limit (and reindroducing the scale a), 
we have 



Suf^ a 2 (d x u x ) 2 
{Au 2 + Auf)^ a 2 Q(0*u x ) 2 + ^{d y u x f 

we find finally our desired result, namely 

H e ff = i j d 2 r [n e s(dyU x ) 2 + K cfi (d x u x ) 2 
with 

„ „ , ; oo L , 9(««) 2 (-, 17 \ 



(B22) 
(B23) 

(B24) 



(B25) 



(B26) 



K 



eff 



where 



K. 



cir 



1 + 



H 2 

64tt 2 



1 -8v- 



23+ 104?; \ k B T 1 
3ko / p 2 Uk J ' 
(B27) 



Meff 



Ka — 1) Voe 



= - (9(ko) 2 + 15/sa + 15) V Q e- Ka 



(B28) 
(B29) 
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